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ABSTRACT 


The  problem  of  electromagnetic  radiation  and  scattering  from  conducting 
bodies  of  revolution,  both  loaded  and  unloaded,  is  considered.  The  solution 
is  obtained  by  the  method  of  moments  applied  to  the  integro-differential 
equation,  and  is  expressed  in  terms  of  generalized  network  parameters. 
Computer  programs  are  given  for  plane- wave  scattering  and  for  radiation 
from  apertures.  The  programs. are  valid  for  solid  bodies,  zero  thickness 
conductors,  and  bodies  with  points  and  edges.  The  solution  for  planar 
bodies  is  also  a  solution  for  diffraction  through  apertures  in  a  conducting 
plane,  according  to  Babinet’s  principle.  Representative  computations  are 
given  for  scattering  and  radiation  from  disks,  washers,  toroids,  hemispheres, 
and  cone- spheres. 
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.  INTRODUCTION  AND  SUMMARY 


The  objectives  of  this  project  were  to  develop  practical  methods  for 
computing  the  generalized  impedance  matrix  applicable  to  the  study  of 
electromagnetic  radiation  and  scattering  from  highly  conducting  bodies  of 
revolution.  Specific  topics  considered  were:  (A)  Development  of  a  com¬ 
puter  program  to  allow  numerical  evaluation  of  the  generalized  impedance 
matrix  for  a  conducting  cone-sphere  of  maximum  dimension  approximately  six 
wavelengths.  (B)  Application  of  the  impedance  matrix  to  determine  the 
scattering  properties  of  a  cone- sphere  with  and  without  loading  elements 
on  the  body.  (C)  Application  of  the  impedance  matrix  to  determine  the 
radiation  and  scattering  properties  of  a  cone-sphere  with  antenna  slots 
on  its  sucface. 

The  project  work  was  performed  during  the  period  1  June  1967  to 
30  June  1969.  The  technical  personnel  working  on  the  project  were  Dr. 

Roger  F.  Harrington  and  Dr.  Joseph  R.  Mautz.  During  the  contract  period 
one  previous  scientific  report  was  written  entitled  "Generalized  Network 
Parameters  for  Bodies  of  Revolution,"  dated  May  1968.^  Also,  one  paper 
entitled  "Radiation  and  Scattering  from  Bodies  of  Revolution"  has  been 

[g] 

published  in  Applied  Scientific  Research.  The  material  in  the  above 
listed  scientific  report  and  paper  also  formed  a  part  of  Dr.  Mautz' s 
Ph.D.  thesis,  Syracuse  University. ^ 

The  following  is  a  summary  of  the  research  performed  during  the  project. 

(A)  A  general  program  was  developed  for  computing  generalized  network 
parameters  for  conducting  bodies  of  revolution.  The  theory  and  original 
program  is  given  in  the  previous  report. ^  An  improved  program  is  given 
in  this  report, 

(3)  Aperture  antennas  in  conducting  bodies  of  revolution  were  analyzed, 
and  computations  made  for  body  currents  and  radiation  patterns.  The  theory 
is  given  in  the  previous  report, ^  and  the  computer  program  is  given  in 
tills  report. 
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(C)  Electromagnetic  scattering  by  conducting  bodies  of  revolution  was 
considered.  The  theory  is  given  in  the  previous  report, ^  ^  and  the  computer 
programs  for  plane-wave  scattering  are  given  in  this  report. 

(D)  Radiation  and  scattering  from  loaded  bodies  of  revolution  was 
studied,  and  computations  made  for  representative  bodies.  The  theory  and 
some  computations  for  cone- spheres  are  given  in  this  report. 

(E)  The  problem  of  determining  the  eigencurrents  for  conducting  bodies 
of  revolution  was  considered,  and  representative  computations  were  made.  De¬ 
tails  are  not  included  in  this  report  because  only  a  preliminary  study  was  made. 

(F)  The  problem  of  synthesizing  and  optimizing  radiation  and  scattering 
patterns  from  loaded  bodies  of  revolution  was  considered,  and  representative 
computations  were  made.  Details  are  not  included  in  this  report  because 
only  a  preliminary  study  was  made. 


II.  METHOD  OF  SOLUTION 

The  problem  is  formulated  in  terms  of  potential  integrals,  and  the 

boundary  conditions  at  the  body  are  applied.  This  results  in  an  integral 

equation,  which  is  then  solved  by  the  method  of  moments.  The  matrices 

obtained  thereby  can  be  expressed  in  terms  of  generalized  network  para- 

[Ul 

meters.  A  general  exposition  of  the  procedure  is  given  in  a  recent 

[5] 

monograph.  The  following  is  a  summary  of  the  theory  as  it  applies  to 
the  present  problem. 

4  sj 

Let  fe  denote  the  known  impressed  field  and  g  the  scattered  field 
due  to  currents  on  the  body.  Then  the  total  field  £  is  the  sum  of  the 
impressed  and  scattered  fields,  that  is. 


£  -  r 


(1) 


The  scattered  field  can  be  expressed  in  terms  of  a  vector  potential  4  and 
scalar  potential  $  as 
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Here  S  is  the  surface  of  the  conductor,  R  is  the  distance  from  a  source 

point  to  the  field  point,  £  is  the  surface  current  on  S,  and  a  is  the  sur- 

m-  ■* - 
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face  charge  on  S.  The  current  and  charge  are  related  by  the  equation  of 

1 

E 

continuity 
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The  boundary  condition  requires  that  the  tangential 

vanish  on  S.  Hence, 

component  of  total  £ 

1 

s 

^tan  ~  ”®tan 

(6) 

where  the  subscript  tan  denotes  tangential  component  on  S.  The  problem 

■ 

can  now  be  stated  succinctly  as 

. 
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where  L  is  the  integro-differential  operator 
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A  solution  of  (7)  gives  the  current  ^  on  3,  Usually  we  are  interested 

in  some  functional  of  £,  which  can  be  computed  once 

^  is  known. 

To  effect  a  solution  by  the  method  of  moments, 

let  the  inner  product 

be  defined  as 
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where  £  and  H  are  tangential  vectors  on  S.  A  set  of  expansion  functions  {jjj} 
is  next  defined,  and  the  current  on  S  approximated  by 

<“» 

3 

where  I,  are  constants  to  be  determined*  Equation  (10)  is  substituted 
J 

into  (7),  which,  because  of  the  linearity  of  L,  reduces  to 

=  <U> 

3 

A  set  of  testing  functions  0^}  is  defined,  and  the  inner  product  of  (11) 
with  each  is  taken.  The  result  is 

Ij  %>  UL>  -  £*>  (12) 

3 

i  =  1,2,3,...  •  The  subscript  tan  has  been  dropped  from  because  the  inner 
product- involves  only  tangential  components.  Wc  now  define  the  generalized 
network  matrices 


CzJ  »  (<  Hj,  ULj  >3  (13) 

CvJ  -  (<  ^  >1  (1*0 

Cl]  -  (ij  (15) 

and  rewrite  the  set  (12)  as 

(El  (I)  =  [V]  (16) 


C 2 ]  is  the  generalised  impedance  matrix,  and  (Y]  =  £ 2 3  “  is  the  generalised 
admittance  matrix.  The  inverse  of  (l6) 


CiJ  -  £y]  (V) 


(17) 
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gives  the  coefficients  1^  of  the  current  expansion  (10),  and  hence  is  an 
approximate  solution  to  the  problem. 

The  impedance  elements  of  (13)  are  explicitly 


^  •  (0«*4j  +  3>j)  ds 


(18) 


The  subscript  0  denotes  that  4.  and  <t> 

J  J 


where  we  have  used  (8)  and  (9), 

are  the  potentials  due  to  £  and  0..  If  the  two-dimensional  divergence 

3  3 

theorem  is  applied  to  the  vector  <t>^  on  a  closed  surface,  the  following 
identity  results; 


St>  •  ds  =  - 


‘  a  <is 
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(19) 


If  y.  is  thought  of  as  a  current,  the  charge  associated  with  it  is 


0 
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*1 


(20) 


Now  (lfi)  can  be  ^r’tten  as 


a  jU 


•  4,  +  0.  $,)  ds 

V  ^  V 


(21) 


Equation  ( 21 )  is  more  convenient  for  computation  than  (iS)  because  the 
gradient  operation  on  $  has  been  eliminated. 


The  above  formulas  have  been  specialised  to  bodies  of  revolution  in 

r  1) 

the  previous  report.  “  Also  given  there  are  representative  computations 
for  radiation  end  scattering  from  unloaded  conducting  bodies  of  revolution. 
Computer  programs  and  further  computations  are  given  in  this  report.  An 
alternative  solution  for  the  case  of  plane-wave  scattering  from  unloaded 
conducting  bodies  of  revolution  has  been  given  by  Andresses . 


The  general  theory  of  loaded  antennas  and  scaitcrers  with 
lumped  loads  is  available  in  the  literature.  "  ' '  An  extension  to  contin¬ 
uously  loaded  surfaces  will  now  be  given.  Lumped  loads  can  be  cnr' * i de red 
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as  a  special  case  of  continuous  loads  localized  to  infinitesimal  sections 
of  the  surface. 

A  continuously  loaded  surface  S  is  defined  as  one  for  which  the  total 
tangential  electric  field  is  related  to  the  current  £  on  S  by  an  impedance 
function  of  position  ^  according  to 

Stan'OC"  ■  <2S> 

Now  is  related  to  the  current  £  by 


4n  *  -L<i> 


(23) 


where  L  is  given  by  (8).  Combining  (72)  and  (23),  and  rearranging,  we  have 


This  reduces  to  (?)  when  ^  «  0,  that  is,  when  S  is  covered  by  a  perfect 
conductor.  We  can  now  reduce  (21*)  to  a  matrix  equation  in  the  usual  way  by 
the  method  of  moments,  obtaining 


(2][ij  »  Cv3  -  ta^Ki] 


(25) 


where  (E),  (v),  tli  are  given  by  (13).  (l*),  (15),  and 

ULJ  *  i%,  ^  (26) 

The  solution  to  (23)  for  the  current  coefficient  matrix  is 


til  *  U  ♦  ZLf4  tv)  (27) 

Side  that  this  solution  is  analogous  to  that  for  two  n-port  networks  con¬ 
nected  ih  scries  with  th*»  voltage  source  ?vj. 

The  impedance  function  X  is  sere  over  these  parts  of  B  covered  by  a 

C 


perfect  electric  conductor.  If  subsectional  expansion  and  testing  functions 
are  used,  many  of  the  elements  of  [z^]  ®ay  bs  when  the  surface  S  is  par¬ 
tially  covered  by  an  electric  conductor.  Ir.  such  cases  the  following  alter¬ 
native  solution  may  be  computationally  faster.  Suppose  [z. ]  has  some  zero 
rows  and  columns.  Let  [z^J  be  the  matrix  obtained  from  [z^]  by  deleting 
all  zero  rows,  [z£]  by  deleting  all  zero  columns,  and  [z£c]  by  deleting 
all  zero  rows  and  columns.  Other  matrices  with  the  same  rows  and/or 
columns  deleted  will  be  identified  by  the  same  superscripts.  Then,  multi¬ 
plying  (25)  by  [Y]  =  [z  X],  and  deleting  the  appropriate  rows  and  columns, 
we  have 

[ir]  =  Cy*  1C v]  -  [Yrcj[zjfcKir]  (28) 

The  solution  of  this  for  [zfc][lr]  is 

Li 

[z£CHlr]  =  Iy*0  +  Y[CrX  [Yr][v]  (29) 

where  Ly£cJ  a  [z£C]  X.  The  solution  is  now  given  by  ( IT)  with  (v]  re¬ 
placed  by 


CV']  *  [V]  ♦  [vLJ 


(yo) 


where  (V^j  is  the  matrix  obtained  by  adding  the  appropriate  seros  to 

(v[]  *  -  U£c]lxrj 

*  -  *  Y^r1  (Yriivj  (>1) 

“The  excitation  is  thus  viewed  as  the  superposition  of  the  i«tp**essed 

voltage  ( V i  plus  the  load  voltage  (v.  J  . 

L 

The  computations  of  the  next  section  for  loaded  antennas  and  scat¬ 
tered  were  made  using  this  second  forau let ion.  A  problem  arises  in  the 
cas<.  of  an  open  circuit,  since  then  e lesson ts  of  may  feecosse  infinite. 

is 
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However,  [Y?c]  is  still  well  defined,  and  may  be  obtained  from  by 

"  Xj 

setting  the  "infinite”  elements  to  some  very  large  number. 


III.  REPRESENTATIVE  COMPUTATIONS 

The  graphs  of  this  section  are  intended  to  illustrate  some  of  the  types 
of  computations  that  we  have  made  during  the  project.  Computer  programs 
for  the  unloaded  antenna  and  scatterer  computations  are  described  in  the 
next  section,  and  listed  in  the  Appendix.  Included  in  the  programs  are 
printer  plot  routines,  so  that  rough  graphs  of  the  computations  are  avail¬ 
able  immediately. 


All  programs  have  been  run  on  spherical  bodies,  and  compared  to  the 
classical  eigenfunction  solutions  for  spheres.  Far-aone  quantities,  such 
as  radiation  fields,  agree  very  closely  (usually  within  a  fraction  of  one 
percent)  with  the  classical  solution.  Iiear-2one  quantities,  such  as  currents 
on  the  body,  also  agree  well  (usually  within  a  few  percent)  with  the  classical 
solution.  Examples  of  such  comparisons  are  given  in  the  previous  report. ^ 

The  programs  have  been  written  to  handle  correctly  the  following  types 
of  bodies:  (A)  solid  conductors,  such  as  spheres,  closed  cylinders,  (B)  sero 
thickness  conductors,  such  as  disks,  open  cylinders,  (C)  bodies  intersect¬ 
ing  the  axis,  such  as  spheres,  disks,  (D)  bodies  not  intersecting  the  axis, 
such  as  toroids,  washers,  (£)  bodies  with  points,  such  as  cone-spheres, 

(F)  bodies  with  edges,  such  as  disks  and  open  cylinders,  and  (G)  bodies  with 
conwrs,  such  as  closed  cylinders.  The  programs  have  not  been  written  to 
handle  multi-body  problems,  but  the  codification  required  is  simple.  The 
accuracy  of  computations  depends  on  the  smoothness  of  thr  body  and  of  the 
excitation.  For  example,  solutions  with  bodies  having  points  or  edges  con¬ 
verge  more  slowly  than  solutions  for  bodies  with  gently  curved  surfaces. 
Solutions  for  localised  aperture  excitation  converge  acre  slowly  than  solu¬ 
tions  for  plane-wave  excitation. 


Solutions  involving  sere- thickness  conductors  in  the  z  *  0  plane  are 
interesting  because  they  arc  also  solutions  to  the  dual  nrohlest  of  apertures 
in  an  infinite  conducting  screen  (Eabinefs  principle)."'  J  in  particular, 
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the  duality  relationship  is 


E 

P 


(32) 


where  J, ,  J  is  the  current  on  the  conducting  plate,  and  E  .  E,  is  the  tan- 
0  P  P  0 

gential  field  in  the  complementary  aperture.  A  number  of  computations  for 

the  circular  disk  and  aperture  were  made,  and  compared  to  available  measure¬ 
ments  .  ^  The  agreement  was  within  the  estimated  experimental  accuracy. 

Figure  1  shows  the  current  on  a  circular  disk  of  diameter  D  =  0.8k,  excited 

by  an  axially  incident  plane  wave,  as  computed  by  the  program  of  Section  V. 

Note  that  J  — >  0  and  J  — >  oo  at  the  edge,  as  theory  predicts.  Figure  2 
P  0 

shows  the  computed  bistatic  radar  cross  section  patterns  for  the  same  disk. 

The  curve  labeled  oX^/k  is  the  E-plane  pattern,  and  that  labeled  /k  is 
the  H-plane  pattern.  (See  the  previous  report^ ^  for  a  discussion  of  the 
notation.)  By  duality,  the  same  patterns  are  also  power  patterns  far  the 
field  diffracted  through  a  circular  aperture  in  an  infinite  conducting  plane 
excited  by  a  plane  wave  at  normal  incidence. 


Figure  3  shows  the  computed  current  on  a  circular  washer  of  outside 
diameter  D  =  2.4k  and  inside  diameter  d  =  0.8k,  excited  by  an  axially  inci¬ 
dent  plane  wave.  Again  J  and  J  have  the  expected  behavior  at  the  edges. 

P  0 

Figure  4  shows  the  computed  bistatic  radar  cross  section  patterns  for  the 
same  disk.  By  duality,  these  patterns  are  also  power  patterns  for  the  field 
diffracted  through  a  circular  annulus  in  an  infinite  conducting  plane  ex¬ 
cited  by  a  plane  wave  at  normal  incidence. 


Figure  5  shows  another  example  of  a  bistatic  radar  cross  section  pattern 
for  a  body  of  revolution  excited  by  an  axially  incident  plane  wave.  It  is 
for  a  toroid  of  mean  radius  R  =  0.54  and  circular  cross  section  of  radius 
r  =  0.12k.  The  previous  report^ ^  showed  bistatic  radar  cross  sections  of 
cone-spheres,  as  well  as  the  currents  on  the  bodies.  These  computations 
are  repeated,  in  the  computer  print-out  of  Appendix  B. 


Figures  6,7;  and  8  show  some  examples  of  radiation  from  rotationally 
symmetric  apertures  in  bodies  of  revolution.  .*■  gure  6  is  the  power  gain 
pattern  for  a  circular  disk  of  diameter  D  =  1.6k  excited  by  a  voltage  across 


Figure  2,  Bistatic  radar  cross  section  for  a  conducting  circular  disk  of  C.S\  diameter  excited 
by  art  axially  Incident  plane  wave.  The  E-plane  pattern  is  labeled  cr*0  and  the  H- 
plane  pattern  a*®. 


excited  by  an  axially  incident  plane  wave. 


Figure  *» .  Bistatic  radar  cross  section  for  a  conducting  circular  washer  of  2-Ux.  outside  diameter 
and  0.6k  inside  diameter,  excited  by  an  axially  incident  plane  wave.  E-plane  and 
K-plane  patterns  are  almost  the  same. 


Figure  5.  Bistatic  radar  cross  section  for  a  conducting  toroid  of  1\  mean  diameter  and  0.24X. 

cross-sectional  diameter,  excited  by  an  axially  incident  plane  wave.  The  E-plane 
pattern  is  labeled  a*6  and  the  H- plane  pattern  0*$ . 


Figure  6.  Power  gain  pattern  for  a  conducting  circular  dick  of  1.6k  diameter,  excited  by  a  voltage 
across  a  narrow  slot  of  0.8\  diameter. 


FEED 


gain  pattern  for  a  conducting  toroid  of  IX  mean  diameter  and  0.2hx  cross-sectional 
,er,  excited  by  a  voltage  across  a  narrow  slot  at  the  maximum  diameter. 


Figure  S.  Power  gain  pattern  for  a  conducting  hemisphere  of  l.k\  diameter,  excited  by  a  voltage 
across  a  snail  annular  slot  at  the  center  of  the  flat  face. 
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a  narrow  slot  at  a  radius  R  =  0.k\.  Figure  7  is  the  power  gain  pattern  for 
a  toroid  of  mean  radius  0.5X  and  cross  sectional  radius  r  =  0.12X,  excited 
by  a  voltage  across  a  narrow  slot  at  the  furthest  point  from  the  z  ;xis. 

Figure  8  is  the  power  gain  pattern  for  a  hemisphere  of  radius  R  =  0.7X, 
excited  by  a  voltage  across  a  small  annular  slot  at  the  center  of  the  flat 
plate.  This  physically,  corresponds  to  excitation  by  an  open-ended  coaxial 
line.  It  is  also  equivalent  to  excitation  by  a  small  z-directed  electric 
dipole  at  the  center  of  the  flat  plate.  Radiation  from  apertures  in  car.p- 
spheres  are  shown  in  the  previous  report. ^  These  computations  are  repeated 
in  the  computer  print-out  of  Appendix  C. 

Figures  9  and  ID  are  examples  of  monostatic  radar  cross  sections  for 
bodies  of  revolution  excited  by  plane  waves  of  arbitrary  direction  of  inci- 
dence.  The  curves  of  c  /X  are  for  ^polarized  incident  waves,  and  the 
curves  of  o^/X2  are  for  polarized  incident  waves.  There  are  no  cross- 
polarized  components  of  backscattered  fields  from  bodies  of  revolution. 

Figure  9  shows  the  monostatic  radar  cross  section  of  a  disk  of  diameter  D  =  0.8\. 
Figure  ID  shows  the  monostatic  radar  cross  section  for  a  hemisphere  of  radius 
R  =  0.7X.  The  computer  print-out  of  Appendix  D  gives  the  monostatic  radar  cross 
section  for  a  cone-sphere, ‘ half  cone  angle  10°,  sphere  radius  0.2X. 

The  computer  print-out  of  Appendix  D  also  shows  the  convergence  of  the 
pattern  as  more  modes  are  included  in  the  solution.  Calling  the  modes  which 
vary  as  the  n-th  modes,  for  accurate  far-field  patterns  one  should  in¬ 

clude  all  modes  for  which 


n  <  C  +1 
-  max 


where  Cmax  =  maximum  body  circumference  in  wavelengths.  A  justification  of 
(33)  can  be  made  in  terms  of  spherical  mode  theory.  Those  modes  having  n 
greater  than  are  cut-off,  that  is,  are  highly  reactive  and  do  not  contri¬ 
bute  materially  to  far- zone  fields. 

Similar  convergence  problems  arise  in  aperture  problems  for  which  the 
excitation  is  not  rotationally  symmetric.  To  illustrate  this,  we  have  com¬ 
puted  the  power  gain  pattern  for  a  cone-sphere  excited  by  a  voltage  across  a 
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usostatlc  radar  cross  section  for  a  conducting  hemisphere  of  l.t»\  diameter,  excited 
'  *  plane  xave  incident  at  an  angle  0.  The  ^polarized  and  0-polarized  patterns  ar* 


slot  at  the  cone-to-sphere  junction,  extending  over  a  90°  angular  sector 
of  the  body.  Figure  11  shows  how  the  power  gain  pattern  converges  as  inodes 
are  added  to  the  solution.  All  curves  are  normalized  to  the  power  radiated 
by  the  entire  aperture,  that  is,  to  the  sum  of  the  modal  powers.  Note  that 
criterion  (33)  is  still  a  good  indication  of  convergence. 

Figure  12  illustrates  a  computation  of  monostatic  radar  cross  section 

vs.  frequency  for  a  conical  sector  of  a  sphere.  Our  computations  (x’s)  are 

[l2] 

compared  to  those  by  Schultz  et_al,  who  used  an  eigenfunction  expansion 
and  nonfinal  determination  of  coefficients.  Also  shown  are  an  approximate 
solution  of  Keller  and  measurements  by  Keys.  The  accuracy  of  our  computa¬ 
tions  is  estimated  to  be  better  tian  can  be  read  from  the  graph. 

Computations  for  the  cone-sphere,  both  as  a  loaded  scatterer  and  as  a 
loaded  radiator,  have  been  made.  The  dimensions  of  the  cone-sphere  and 
points  of  loading  are  pictured  in  Fig.  13.  Computations  were  made  for  four 
types  of  loads  as  follows: 

(A)  Short  circuit 

(3)  Open  circuit 

(C)  Resonant  load 

(D)  {•latched  load 

For  the  scattering  problem,  these  locds  were  applied  one  at  a  time  to  each 
of  the  four  loading  points.  The  plane-wave  illumination  was  axially  incident 
on  either  the  point  or  the  sphere.  For  the  radiation  problem,  the  excitation 
was  taken  at  one  of  the  four  points,  and  one  of  the  other  points  was  loaded 
by  one  of  the  above  loads.  We  have  a  complete  set  of  these  patterns  for  all 
permutations  of  the  excitations  and  loads.  Representative  acts  of  scattering 
and  radiation  patterns  are  shown  in  Figs.  l«*  and  13  to  illustrate  the  results. 


iv.  tngaguujgsp  jmtwonK  matrices 


This  program  computes  the  Csl  and  It]  matrices  for  bodies  of  revolution. 
Hie  theory  is  that  of  the  previous  report' *  modified  to  allow  unequal  seg¬ 
ments  of  contour.  The  input  data  consists  primarily  of  Jip  points  specifying 
the  contour  plus  the  afiede  number-  The  output  consists  primarily  of  the 


Figure  12.  Monostatic  radar  cross  section  for  a  conical  sector  of 
a  sphere,  as  computed  and  meas'  red  by  various  persons . 
Our  computations  are  shown  by  x's. 


Figure  13.  Conducting  cone-3phere  and  points  of  loading 


INCIDENT 


Figure  1*.  Bistatic  scattering  patterns  (c/\  vs.  e)  for  a  Loaded  cone-sphere,  Loads  are  resonant 
bands  at  (a),  (b),  (c),  or  (d).  The  plane  wave  is  axially  incident  on  the  tip,  and  the, 
patterns  are  in  the  E  plane.  Cone  angle  is  20  ,  sphere  diameter  is  0.4\.  (Maxiirna  o/V 
for  unloaded  cone-sphere  is  O.38. ) 


Figure  15.  Power  gain  patterns  (G  vs.  q)  for  a  loaded  cone-sphere  slot  radiator.  Feed  is  at  (b). 

loads  are  resonant  bands  at  (a),  (b),  (c),  or  (d).  Cone  angle  is  20°,  sphere  diameter 
is  O.U\.  There  is  no  0  variation. 
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elements  of  the  generalized  admittance  matrix  for  the  particular  mode.  The  pro¬ 
gram  listing  plus  sample  input-output  data  are  given  in  Appendix  A. 

The  subroutine  LINEQ(LL,C)  replaces  the  LL  by  LL  matrix  C  by  its  inverse.  C 

till 

is  a  complex  one-dimensional  variable.  The  (l,j)  matrix  element  resides  in 
C(l  +  (J-1)  *  LL).  If  LL  >  58,  the  dimension  of  the  variable  LR  must  be  increased 
to  at  least  LL. 

Punched  card  data  is  read  early  in  the  main  program  according  to 

50  READ(1.51,  END  =  52)  NN,  NP,  NPHI,  BK 

51  FOmT(313,  Ell. 7) 

READ(l,53)(RH(l),  I  «  1,  HP) 

READ(l,53)(ZH(l),  I  =  1,  NP) 

53  FORmT(lOF8.1) 

The  excitation  will  haye  dependence,  NPHI  has  the  same  meaning  as  on  page  55  of 

the  previous  report. ^  BK  is  the  propagation  constant  k  =  .  An  odd  number 

NP  of  points  are  taken  from  the  generating  curve  of  the  body  of  revolution.  KH(l) 
and  ZH(l)  are  respectively  the  distance  p  from  the  axis  (z  axis)  of  the  body  of 
revolution  and  the  corresponding  z  coordinate.  RH(l)  may  be  zero  only  when  1=1 
or  I  =  NP.  If  the  generating  curve  closes  upon  itself,  care  must  be  taken  to  make 
the  coordinates  at  I  =  1  identical  to  those  at  I  =  NP.  The  data  is  printed  just 
after  it  is  read. 

If  the  first  and  NP^h  data  points  are  the  same,  the  second  and  third  data 
points  are  overlapped  into  the  NP  +  1  and  NP  +  2  position  and  NP  is  increased  by 
two  before  going  to  do  loop  57  which  interpolates  to  find  the  distance  DH(l)  be¬ 
tween  the  I  ^  and  (i  +  l)*'*1  data  points,  the  radius  p  midway  between  the  1^  and 
t/ll 

(I  +  1)  data  points,  the  cori'esponding  z  coordinate  ZS(l),  and  the  sine  SV(l)  and 

til 

cosine  CV(l)  of  the  angle  between  the  z  axis  and  the  straight  line  from  the  I  to 

the  (I  +  1)  data  points.  Next,  a  few  constants  are  entered  among  which  PI  is  n 

and  ETA  is  the  intrinsic  impedance  tj  =  376.707  taken  from  the  second  edition  of 

Smythe.^^  The  third  edition  of  Smythe  gives  -q  =  376.730.  DO  loop  117  computes 

the  t  coordinate  (arc  length  along  the  generating  curve)  TJ(I)  of  the  peak  of  the 
th 

I  triangular  expansion  function.  TJ  is  printed.  DO  loop  K>  calculates 
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CSM(K+(l-NN+l)*NPHl)  =  *  COS 

and  I  =  NN-1,  NN,  NN+1. 

DO  loop  1 6  calculates 


mi 


for  K  =  1,2,. ..NPHI 


dh'(j) 


/  dt 


/  d 0 
0 


3-jkR 

"kR 


cos  n 0 


(34) 


The  t  integration  is  over  the  portion  (assumed  to  be  a  straight  line)  of 

the  generating  curve  between  the  j  and  (J  +  l)th  data  points,  0  is  the 

angle  swept  when  the  generating  curve  is  revolved  about  the  z  axis.  R  is 

the  distance  between  some  field  coordinate  (t',  0'  =  0)  and  the  source 

coordinates  (t,  0).  The  expression  (34-)  is  proportional  to  the  scalar 

potential  at  (t',  0)  coming  from  a  ^  charge  density  on  the  lateral 

P 

surface  of  the  frustum  of  a  cone.  The  index  J  of  do  loop  1 6  indicates  that 

the  source  coordinate  t  is  between  the  <j  and  (j  +  l)^n  data  points.  The 

index  I  of  do  loop  17  indicates  that  the  field  coordinate  t*  is  midway  be- 

tween  the  I  and  (i  +  1)  data  points.  Except  when  the  field  point  lies 

on  the  surface  charge,  the  surface  charge  is  approximated  by  line  charges 

at  0  =  >  K  =  1,2,... NPHI.  When  the  field  point  lies  on  the  surface 

charge,  the  surface  charge  -cos-n&  is  approximated  by  cos  ( .5*Pl/NPHl)/R(l) 

PI  p  jt 

when  0  <  0  <  ,  but  for  the  region  <  0  <  x  the  line  charges  at 

( k-  .5  )*Pl NPHI 
0  =  —  ,  K  =  2,3,  ...NPHI  are  used.  Do  loop  5  stores  the  results  of 

wrni  (k-.5^*PI 

the  t  integrations  at  0  =  ,  K  =  1,2,  •••NPHI  in  GS,  Do  loop  13 

multiplies  the  GS(K)  by  j— j  *  cos  ~)  and-  sums  them  over  K  to 

obtain  expression  (34)  stored  in  G.  Do  loop  68  obtains  n  =  NN-1,  NN,  and 
NN+1.  The  location  G(l  +  (j-l)*(NP-.l)  +  (NP-1)  *  (NP-l)  *(n-NN+l))  cor¬ 
responds  to  field  point  I,  source  region  J  and  mode  n.  If  the  generating 
curve  closes  upon  itself,  the  present  NP  is  two  more  than  the  original  NP 
of  the  input  data. 

The  impedance  matrix  2  will  be  the  same  as  that  of  equation  (45)  of 
tue  previous  report^  except  that  the  previously  defined  G^'s  will  be  used 


and  the  T  ani  T'  are  modified. 
P  P 


Also,  since  the  new  Gn's  are  divided  by  k 


to  make  them  insensitive  to  the  absolute  si2e  of  the  body  of  revolution,  the 
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new  Z's  must  be  multiplied  by  k.  Because  of  possible  unequal  spacing  of 
data  points  the  8  numbers  T^  and  were  stored  separately  for  each  tri¬ 
angular  'expansion  function.  In  do  loop  74,  the  index  J  refers  to  the  Jth 

triangular  function  while  the  numbers  T  and  Tf  are  stored  in  T  and  TP 

th  P  P 

respectively.  The  J  triangular  function  is  approximated  by  4  pulses  of 
base  lengths  DH(l  +  2  *(J-1)),  I  =  1,2, 3,4.  The  derivative  of  the  Jth  tri¬ 
angular  function  is  Just  a  pulse  doublet.  Denoting  base  lengths  by  A^, 
by  aM 


T(1  +  4*(j-l))  = 
T(2  +  4*(j-l))  = 
TP(1  +  4*(j-l))  = 
TP(2  +  4*(j-l))  , 


2*(A1+^) 

(a1  +  .5*^)*^ 

+  ^  - 

^  "■ 

*l  +  ^> 


(35) 


Similar  expressions  involving  A^  and  A^  are  used  to  treat  the  trailing  edge 
of  the  triangular  function.  Expressions  (35)  are  functional  values  (T  for 
the  triangle  and  TP  for  the  derivative  of  the  triangle)  weighted  by  the  base 

length  in  question.  The  base  length  weighting  is  necessary  because  the  G  's 

11  n 
have  a  factor  =  y  which  must  be  offset.  Also,  the  factor  A.  is  neces¬ 

sary  for  the  field  integration.  For  the  field  integral,  the  triangulaz' 
function  and  its  derivative  are  each  approximated  by  four  delta  functions. 

Unless  the  generating  curve  closes  upon  itself,  there  will  be  no  peak 
of  triangular  function  at  the  starting  point  and  end  point.  The  expansion 
function  which  is  actually  a  triangular  function  over  radius  p  will  reduce 
to  a  constant  at  a  starting  point  or  end  point  on  the  z  axis.  However,  if 
the  starting  point  or  end  point  is  not  on  the  z  axis,  the  expansion  functions 
reduce  to  aero  there,  in  problems  involving,  perhaps,  a  conducting  disk,  the 
electric  current  normal  to  the  rim  goes  to  zero  as  the  square  root  of  the 
distance  from  the  rim  and  the  current  parallel  to  the  rim  becomes  large  as  one 


over  the  square  root  of  the  distance  from  the  rim,  at  least  when  an  incident 
plane  wave  propagates  perpendicular  to  the  plane  of  the  disk.  The  triangular 
functions  used  to  expand  the  0  directed  current  were  modified  near  the  ex¬ 
tremities  of  the  generating  curve  where  no  overlapping  occurs  remote  from 
the  z  axis.  The  triangular  function  (associated  with  the  0  directed  current) 
nearest  such  an  extremity  was  modified  by  increasing  its  value  at  the 
extremity  from  zero  to  two.  The  logic  from  statements 
115  to  ll6  performs  this  modification.  Whereas  the  triangular  functions 
associated  with  the  t  directed  current  were  stored  in  T,  the  functions 
associated  with  the  0  directed  current  are  assigned  TR. 

Using  Gn,  T,  TP,  and  TR  the  impedance  matrix  is  calculated  inside  do 
loop  30.  The  calculation  of  the  present  impedance  elements  is  very  similar 
to  a  previous  calculation  using  equation  (L5)  of  the  previous  report. 

Statement  8l  CALL  LINEQ  (NM2,  Z)  inverts  the  impedance  matrix  stored 
in  Z  to  obtain  the  admittance  matrix.  The  admittance  matrix  is  printed  in¬ 
side  do  loop  93.  If  superscripts  t0  are  replaced  by  rs  and  if  r  =  0  denotes 
t,  r  =  1  denotes  0  and  similarly  for  s,  then  (Zn  )„,  i  =  1,2,  •••N-l  is 

given  by  the  (j  +  r  *  (N-l)  +  s  *  2  *  (N-l)]th  execution  of  write 

statement  96.  Here,  N-l  is  the  number  (either  or  ^2~)  of  triangular 
functions.  Immediately  following  statement  8l,  the  admittance  matrix  of  the 
array  Z  is  recorded  on  a  direct  access  device  by 

WRITE (6)  (Z(I),  I  =  1,  NZ) 

If  no  direct  access  device  is  available,  statements  should  be  added  to  either 
punch  the  admittance  matrix  on  cards  or  store  it  on  tape.  At  the  end  of  the 
program,  execution  is  transferred  back  to  read  statement  5C.  If  there  is 
another  set  of  data,  it  is  processed.  If  not,  execution  terminates. 

The  dimension  statements  in  the  main  program  will  not  accommodate  NPHI 

larger  than  to  nor  NP  larger  than  t 1  if  the  generating  curve  closes  on  itself. 

If  the  generating  curve  docs  not  close  on  itself,  NP  may  be  as  large  as  ^3. 
Assuming  that  the  generating  curve  might  close  on  itself,  minimum  dimensions 


COMPLEX  Z((NP-l)*(NP-l)),  GS(NFHI),  G(2*(NP+1)*(NF+1) ) 
DIMENSION  RH(NP+2),  ZH(NP+2),  DH(NP+2), 

DIMENSION  SV(NPH),  CV(NP+1),  ZS(NP+1),  H(NPH), 

ANG(NPHI),  AC(NFHl),  CSM(3*NPHI) 

DIMENSION  TP(2*(NP-l)),  T(2*(NP-L)),  TR(2*(NP-l)) 

The  NP  referred  to  above  is  the  original  NP  of  the  input  data. 


.  PLANE-WAVE  SCATTERING.  AXIAL  INCIDENCE 


This  program  computes  the  current  on  a  conducting  body  of  revolution 
and  scattering  patterns  for  excitation  by  a  plane  wave  axially  incident  on 
the  body.  The  input  consists  primarily  of  NP  points  defining  the  body  con¬ 
tour,  plus  the  [y]  matrix  for  the  n  =  1  mode.  The  output  consists  primarily 
of  the  current  on  the  body,'  the  scattered  far- zone  field,  and  the  bistatic 
radar  cross  section.  Printer  plot  routines  are  included  to  graph  both  the 
current  and  the  radar  cross  sections.  The  program  listing  plus  sample  inpi^t- 
output  data  are  given  in  Appendix  B. 


Subroutine  PLANE(WR,  THR,  NT)  provides  the  measurement  matrix  elements 
of  equations  (77)  and  (8l)  of  the  report. ^  NT  angles  (in  radians)  appear 
as  input  in  THR.  Other  input  appears  in  the  common  statement.  U  is  the  com¬ 
plex  constant  (0.,1. ),  R,  ZS,  SV  and  CV  are  respectively  the  radius  p,  the  z 
coordinate,  and  the  sine  and  cosine  of  the  angle  between  the  direction  of 
the  generating  curve  and  the  z  axis.  R(l),  ZS(T),  SV(l)  and  CV(l)  are 
evaluated  midway  between  the  and  (i+l)^1  data  points  RH  and  ZH  by  draw¬ 
ing  a  straight  line  between  them.  BK  is  the  propagation  constant,  there  are 
(NP-l)  R’s,  the  excitation  has  e^^^  dependence,  and  T  and  TR  are  the  pre¬ 
viously  computed  weighted  samples  of  the  triangular  functions  for  both  the 


t  and  <f>  directed  electric  currents.  When  execution  returns  to  the  main  pro¬ 
gram,  the  (Rn)i  of  equation  (77)  or  (8l)  of  the  report^  will  be  stored  in 
WR(i  +  (K~1.)*NM  *  (L-l)*t*NM)  where  K  =  1,2, 5,h  denotes  respectively 


<Ol’ 

the  number 


(R^).  and  (rJ*)  .  As  defined  in  the  subprogram,  NM  is 

of  the  triangular  function's  and  L  denotes  the  Lwl  of  the 


angles  e  . 
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In  the  subroutine  PLANE,  do  loop  153  puts  J!  in  FK(j+l).  Do  loop  156 
sweeps  the  NT  angles  e  >  Do  loop  302  computes  the  required  Bessel  functions 

y  (-if  (|)n+2K 

according  to  Jn(x)  =  Z^_ - KI  (n+K) !  *  If  the  series  for  Jn(x)  does  not 

converge  readily,  do  loop  155  may  be  satisfied  in  which  case  error  stop  155 

is  reached.  When  p  =  R(j),  (k  sin  e^)  will  be  stored  in 

Bj(j  +  (K+l)*(NP-l) )  for  K  =  -1,0,1.  To  avoid  direct  computation  of  a 

Bessel  function  of  negative  order  when  NN  =  0,  J  ^(x)  is  obtained  from  J^(x) 

in  do  loop  309*  With  pf.  being  represented  by  either  T  or  TR,  the  t  inte- 

1  r  ii 

gration  of  equations  (77)  and  (8l)  of  the  report  is  simulated  in  do  loop  301. 
The  pf.  triangular  function  is  approximated  by  four  delta  functions. 

^  V 

If  NP  (either  the  original  NP  of  the  data  or  2  plus  it  if  the  generat¬ 
ing  curve  overlaps)  is  larger  than  43,  the  dimension  of  some  of  the  variables 
in  PLANE  must  be  increased.  These  variables  and  their  dimensions  are: 

C0MM3N  R(NP-l),  ZS(NP-l),  SV(NP-l),  CV(NP-l), 

T(2*(NP-3 ) ), TR(2*(NP-3 ) ) 

DIMENSION  BJ(3*(NP-1)) 

The  common  statement  in  the  subroutine  PLANE  should  be  identical  to  that  in 
the  main  program.  In  the  unlikely  case  that  the  mode  number  NN  is  larger 
than  17,  the  dimension  of  FK  must  be  increased  to  NN+3.  Since  the  arguments 
WR  and  THR  are  only  dummy  variables,  their  present  dimension  of  1  is  always 
sufficient. 

Tlie  subroutine  REQRD(K1,  K3,  L)  rearranges  the  first  L  elements  of  K3 
in  descending  order.  Also,  the  first  L  elements  of  Kl  are  rearranged  so 
as  to  maintain  the  original  correspondence  between  the  elements  of  K3  and 
those  of  Kl. 

Punched  card  data  is  read  in  early  in  the  main  program  according  to 

50  READ(1,51,  END  =  52)  NN,  NP,  NT,  BK 

51  FORMAT(3I3,  El1*.?) 

RKAD(L,53HRH(I),  I  »  1,NP) 

READ(1,53)(ZH(I),  I  »  1,NP) 

53  F0RMAT( 10F8.4 ) 
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All  the  input  variables  except  NT  are  the  same  as  the  variables  of  the 
previous  program.  The  receiver  angles  $ are  given  by 


rt(i-l) 

NT-1 


i  =  1,2, ♦••NT 


The  punched  card  data  is  printed  out  immediately  after  it  is  read  in. 

Statement  126  checks  to  see  whether  the  first  and  last  points  on  •the 
generating  curve  are  identical.  If  they  are  identical,  the  variable  KL  is 
set  equal  to  0  and  RH  and  ZH  are  expanded  by  overlapping  the  second  and 
third  points  on  the  generating  curve.  In  this  way,  an  expansion  function 
centered  about  the  last  original  data  point  is  obtained.  If  the  first  and 
last  original  data  points  are  not  identical,  KL  is  1  and  execution  proceeds 
to  statement  58  in  which  case  there  are  no  expansion  functions  centered 
about  the  first  and  last  data  points.  The  logic  between  statements  58  and 
78  prepares  input  data  for  the  subroutine  PLANE.  Statement  85  inserts  the 
receiver  voltage  matrices  into  WR.  Statement  127  inputs  the  admittance 
matrix  appearing  in  (17)  of  the  report. ^  This  admittance  matrix  will  have 
been  generated  by  the  previous  computer  program. 


Nearly  all  the  statements  following  statement  127  are  enclosed  in  do 

loop  108.  When  INC  =  1,  the  electric  field  of  the  incident  plane  wave  is 

Wher.  inc  <=  2  it  is  -i^e^2.  The  I  column  vector  of  equation 

(17)  of  the  report  is  computed  in  do  loop  2  and  stored  in  TI.  The  inner 

do  loop  3  uses  transmitter  voltage  matrices  obtained  from  the  last 

(0^=  jt,  INC  1)  and  first  ^0  =0,  INC  =  2)  receiver  voltage  matrices. 

Equation  .(82)  of  the  report'1^  and  the  fact  that  is  even  in  n  and  R?^ 

is  odd  in  n  are  used  to  relate  the  transmitter  voltage  matrices  to  the 

receiver  voltage  matrices.  Do  loops  7  and  8  compute  TJ(J),  the  arc  length 

till 

which  locates  the  center  point  of  the  J  expansion  function  along  the 
generating  curve.  The  total  length  of  the  generating  curve  is  normalised 
to  10.  Do  loop  9  performs  the  matrix  multiplication  indicated  in  equation 
(65)  of  the  report.1  i  The  measurement  is  stored  in  El  and  E2  correspond¬ 
ing  to  the  b  and  $  receiver  polarizations.  In  do  loop  ii,  El  and  E2  are 

2 

normalised  so  that  the  quantity  o/\  associated  with  each  one  is  the  square 
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of  its  magnitude.  The  phase  of  El  is  set  equal  to  zero  in  the  backscatter- 

ing  direction  which  is  =  n  for  the  e  ^ z  incidence  and  6r  =  0  for  the 
e+jkz  incj_^ence< 

Do  loop  128  divides  the  expansion  function  coefficients  by  p  to  obcain 
the  current.  A  multiplication  by  2  or  2j  converts  from  the  exponential  e*^ 
to  cos  0  or  sin  0.  The  factor  tj  is  necessary  to  normalize  the  current  to 
the  incident  magnetic  intensity  g1.  Do  loop  129  averages  the  0  directed 
current  J  according  to 


i-1)  4-  2J^(i)  +  J^H-l) 


The  averaging  was  deemed  necessary  to  attenuate  extraneous  oscillations  in 
J^.  If  the  generating  curve  does  not  close  upon  itself,  the  first  and  last 
are  not  averaged.  Do  loop  4  prints  the  arc  length  and  real  and  imaginary 
parts  and  magnitude  of  the  t  and  0  directed  current.  For  an  incident  mag¬ 
netic  intensity  g1  =  (INC  =  1)  or  H1  =  lys-3*2  (INC  =  2),  the  current 

sinor  H  is  given  in  terms  of  the  printed  output  JT  and  J 0  by 


s I  =  n^JT  cos  0  J0  sin  A 

Do  loop  12  print s  the  scattering  patterns  as  functions  of  0  .  SIG  X0  is  o/i? 

O  1'  _ 

in  the  E  plane  and  SIG  X0  is  o/K  in  the  H  plane.  LSIG  X0  is  log.  o/k2 

O 

in  the  E  plane  and  DSIG  X0  is  o/\  in  the  H  plane.  SIG  X0  and  SIG  X0 

are  respectively  the  magnitudes  squared  of  SX@  and  SX0 .  The  phase  of  SXq 
is  taken  to  be  zero  in  the  bockscattering  direction.  The  magnitudes  (MAG) 
and  phase  (AUG)  of  SX0  and  SX0  are  also  tabulated.  The  scattered  elective 
field  will  be  proportional  to 


\iQ  (sxe)  COG0  *■  (3X0)  S1.50 


All  the  statements  following  do  loop  12  are  devoted  to  obtainixrg  plots 
of  the  previously  tabulated  quantities  MAG  JT,  MAG  J0,  SIG  X§,  and  SIG  X0. 
Do  loop  12  scales  the  current  for  plotting.  K3  is  the  t  dii-ectcd  current 


and  K4  is  the  0  directed  current.  K3  has  been  assigned  the  abscissa  K1  and 
Kl  the  abscissa  K2.  After  execution  of  statements  ll  and  15,  K3(<J)  will 
be  the  <j  n  largest  K5  and  Kl(j)  its  abscissa.  The  quantities  K5  and  Kl 
with  abscissas  K1  and  K2  are  plotted  by  do  loop  20.  Entering  do  loop  20, 

K5  =  1*  One  line  is  printed  for  each  J.  If  Kj(K5)  is  greater  than  or 
equal  to  (51-J),  an  X  is  written  after  Kl(K5)  blanks  and  K5  is  incremented 
by  one.  If  K3(K5)  is  still  greater  than  or  equal  to  (51-J),  another  X  is 
written  after  Kl(K5)  blanks  and  K5  is  incremented.  The  process  continues 
until  K3(K5)  is  less  than  (51-J)  in  which  case  statement  20  is  reached. 
Similar  action  is  associated  with  Kk  and  K2.  If  K3(NM)  is  larger  than  or 
equal  to  zero,  K3(NM+1)  will  be  observed,  but  K3(NM+1)  has  been  set  equal  to 
-1  immediately  after  do  loop  ll. 

The  patterns  SIG  X9  and  SIG  X0  are  plotted  in  a  similar  fashion. 
Actually,  the  logarithms  of  SIG  XS  and  SIG  X0  are  plotted  but  the  scale  is 
graduated  in  terms  of  SIG  X0  and  SIG  X0  to  simulate  semi-log  paper.  Do  loop 
80  generates  the  abscissas  K1  and  K2  and  ordinates  K3  and  K4.  Statements 
l6  and  I?  arrange  K3  and  Kl  in  descending  order.  Do  loop  87  plots  SIG  X9 
and  SIG  X0  as  X's  and  0's  respectively.  After  do  loop  108  is  satisfied, 
execution  goes  back  to  the  read  statement  50.  If  another  set  of  input  data 
is  encountered,  execution  traverses  the  program  again.  If  not,  execution 
stops. 

Let  NP  be  either  the  original  NP  of  the  data  if  the  generating  curve 
does  not  close  upon  itself  or  2  plus  the  original  NP  of  the  data  if  the 
generating  curve  closes  upon  itself.  When  NP  >  13  or  NT  >  75,  some  variables 
require  larger  dimensions.  These  variables  are  listed  along  with  their 
minimum  dimensions. 


COMPLEX  Y(l*NM«!iM),  WR(1*NT»NM),  TI(2*8M), 

L‘3(NK),  Ki(NT),  E2(NT) 

CObSm  R(NT-l),  XS(NP-l),  SV(SP-l),  CV(NP-l),  l(l*NK),  TR(1*NM) 

DIMENSION  RH(liP),  ZK(NP),  DH(NP-l),  TJ(*«K),  THR(ltT) 

DIVISION  Kl(NM),  XI (NT),  K2(NK),  K2(NT),  K3(N«),  K3(ST),  Kl(NM),  Kl(ST) 
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The  variables  Kl,  K 2,  K3  and  K4  must  be  dimensioned  at  least  as  large  as 

NP-3 

the  larger  of  NM  and  NT.  Here,  as  in  the  main  program,  NM  =  — ^  . 

VI.  APERTURE  RADIATION.  ROTATIONAL  SYMMETRY 

This  program  computes  the  current  on  a  conducting  body  of  revolution 
excited  by  one  or  more  rotationally  symmetric  apertures.  It  is  assumed 
that  the  tangential  electric  field  is  known  over  the  apertures.  The  input 
consists  primarily  of  NP  points  defining  the  body  contour,  the  [Y]  matrix 
for  the  n  =  0  mode,  plus  tangential  E  at  the  NP  points  on  S.  The  output 
consists  primarily  of  the  current  on  the  body,  the  radiation  field,  and 
the  normalized  power  gain  patterns.  Printer  plot  routines  are  included  to  graph 
both  the  current  and  the  power  gain  patterns.  The  program  listing  plus 
sample  input-output  data  are  given  in  Appendix  C. 

The  program  is  similar  to  the  one  concerning  scattering  for  axially 
incident  plane  waves  except  that  the  excitation  voltage  matrix  is  deter¬ 
mined  by  input  data  instead  of  by  an  incident  plane  wave.  Also,  a  gain 
instead  of  a  scattering  cross  section  is  obtained  from  the  measurement. 
(Equation  (69)  of  the  report. 

The  subroutine  PLANE  provides  the  measurement  matrix  elements  of 

f '  1 

equations  (??)  and  (Si )  of  the  report1 for  the  special  case  n  =  0.  This 
subroutine  PLANE  is  similar  to  the  subroutine  of  the  same  name  used  in  the 
program  concerning  scattering  for  axially  incident  plane  waves.  For  the  case 
n  0,  and  may  be  dispensed  with  since  they  are  zero. 

The  subroutine  REOHD  is  exactly  the  same  as  the  one  of  the  same  name 
appearing  in  the  program  concerning  scattering  for  axially  incident  plane 
waves . 

Punched  card  data  is  read  early  in  the  main  program  according  to 
READ ( 1,51, END  *  52)  KK,  NP,  NT,  BK 

R£AS(i,53)(RK(l),  I  »  1,  NP) 

R£AD(l,53)(23i(l),  1  *  1,  NP) 
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53  format(iof8.4) 

IF(KX.EQ.2)  GO  TO  40 
KEAD(l,68)(El(l),  I  =  1,NP) 
168  FORMA.T(7E11.4) 

IF(KK.Eft.l'  GO  TO  41 

40  READ(l,l68)(E2(l),I  =  1,  HP) 

41  KL  =  1 


KK  =  1  if  the  aperture  electric  field  has  only  a  t  component. 

KK  =  2  if  the  aperture  electric  field  has  only  a  <p  component. 

KK  =  5  if  the  aperture  electric  field  has  both  t  and  <p  components.  NP  points 

define  the  generating  curve  of  the  body  of  revolution.  The  first  point  and 
the  point  mark  the  extremities  of  the  generating  curve.  NP  must  be 

odd.  The  gain  pattern  will  be  observed  at  NT  angles  0^.  In  degrees. 


& 

r 


1,2, •••  NT. 


BK  is  the  propagation  constant  k  3  w/yT  •  Rli( X )  and  2H(l)  are  respectively 
the  cylindrical  coordinate  radius  p  and  axial  distance  z  at  data  point  I  of 
the  generating  curve.  £1  and  £2  are  complex.  If  KK  4  2,  El(l)  is  the  t 
component  of  the  axially  symmetric  aperture  electric  field  at  point  I  on  the  gen¬ 
erating  curve.  All  the  Si's  except  Kl(J)  being  cero  represents  an 
aperture  field  between  the  ( J- 1 ) and  (j+i)^1  points  or  between  the 

f  h  **  ii,  tj, 

J "  and  (d*i)  *  points  if  J  »  1  or  between  the  (J-i)  and  J  points  if 
J  *  NP.  If  KK(1)  *  RH(SP)  4  0  and  2K ( 1 }  ■  3H(NP},  an  Ki(  1}  aust  be  ac¬ 
companied  by  an  identical  Sl(NP).  If  KK  4  1,  £2(1}  represents  a  0  directed 
aperture  field  in  exactly  the  same  way  as  Sl(l}  represented  a  t  directed 
aperture  field.  If  KK  «  3,  both  SI  and  S2  occur  simultaneously.  The 
punched  card  data  is  printed  out  aeon  after  it  is  read. 


Statement  126  sets  KL  *  Q,  overlaps  two  of  the  data  points  arid  adds  2 
to  UP  if  the  generating  curve  closes  upon  itself.  Do  loop  5?  draws  a 
straight  line  between  the  (!♦!)  *  and  I  data  points  to  obtain  J3K(l},  R(l),  2o(l), 


S V  ( I ) ,  and  CV(l),  the  increment  of  arc  length,  the  cylindrical  coordinate 
radius,  the  axial  distance,  and  the  sine  and  cosine  of  the  angle  between 
t  and  ,he  z  axis,  all  evaluated  midway  between  the  (i+l)^  and  1^  data 
points.  Do  loop  1  stores  the  NT  angles  0^  in  THE. 

The  triangle  functions  (equation  31  of  the  report^ )  are  modified 
slightly  if  the  data  points  on  the  generating  curve  are  not  exactly 
equally  spaced. 


T.(t) 


where  i  = 


Ltll 
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til 

The  subscript  j  on  t  indicates  the  J  data  point.  For  J  odd,  this 

*r  4* 

(-T-)  1  function  is  characterized  by  4  numbers  T(2*J-6+l),  I  =  1*2,3,'*. 

_  y_ 

Coming  out  of  do  loop  T(2*J-6^l)  is  T.(t)  midway  between  the  (j-3+,I) 

4-h 

and  (<J-2+I)'J  data  points  multiplied  by  the  arc  length  included  between  the 
(j-3+l)ti:i  and  (J-2+l)tl1  data  points.  The  multiplication  by  the  increment 
of  arc  length  factor  facilitates  integrations  over  the  source  and  field 
regions . 
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If  KL  4-  0  and  if  RH(l)  is  not  zero,  the  body  of  revolution  has  a  rim. 

The  previously  defined  T(t)  nearest  the  rim  will  be  zero  at  the  rim.  Since 
the  0  directed  current  is  not  expected  to  be  zero  at  the  rim,  the  T(t) 
associated  with  the  0  directed  current  nearest  the  rim  was  changed  to  2  at 
the  rim.  This  amplitude  2  makes  the  function  linear  in  t  when  the  points 
tj  in  question  are  equally  spaced.  The  logic  between  statements  74  and  j8 
represents  the  functions  associated  with  the  0  directed  current  by  TR. 

'  Do  loop  7  computes  the  arc  length  TJ  at  every  other  data  point  along 
the  generating  curve.  Do  loop  8  normalizes  TJ  for  a  generating  curve  of 
length  10. 

Using  T  and  TR  to  invoke  a  pulse  approximation  to  the  T(t)  functions, 
do  loop  44  uses  the  aperture  electric  field  data  to  calculate  and  store  in 
E3  the  excitation  voltage  matrices,  (equation  (27)  of  the  report^).  If  E3(j) 
represents  excitation  in  the  t  direction,  E3(J+NM)  represents  a  correspond¬ 
ing  excitation  in  the  0  direction.  Actually  if  KX  4  5.  there  is  only  one 
type  of  excitation  and  either  E3(j)  or  E3(-T+NM)  is  automatically  set  equal 
to  zero. 

Statement  183  reads  the  admittance  matrix  Y  from,  in  this  case,  direct 
access  data  set  6.  Do  loop  134  calculates  the  coefficients  TI  of  the  current 
expansion  functions  by  multiplying  the  admittance  matrice  Y  by  the  excitation 
voltage  matrix  E3.  If  Tl(j)  represents  a  t  directed  current,  then  Tl(j+NM) 
represents  the  corres ponding  0  directed  current.  Since  the  admittance  sub¬ 
matrices  Y^  and  Y^  are  zero,  the  t  directed  aperture  field  will  excite 
only  a  t  directed  electric  current  and  the  0  directed  aperture  field  only 

4 

a  0  directed  current .  Ftili  in  do  loop  Ip4,  the  power  required  to  main¬ 
tain  the  aperture  field  is  obtained  by  summing  the  products  of  the  complex 
conjugate  of  a  coefficient  of  a  current  expansion  function  with  the  corres¬ 
ponding  excitation  voltage  matrix  element.  The  power  associated  with  the  t 
directed  aperture  electric  field  is  put  in  PI  and  the  power  associated  with 
the  0  directed  aperture  electric  field  In  P2.  Since  only  the  real  part  of 
the  complex  power  is  required  to  normalize  the  gain,  neither  Pi  nor  P2  are 
defined  to  be  complex  variables.  PI  and  F2  are  written  just  after  do  loop 
104.  The  sample  printer  output  for  the  cone  sphere  data  came  from  an 
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earlier  version  of  the.  program  which  printed  PI  and  P2  according  to 

WRITE(3,l69)  PI,  P2,  P3 
169  F0RMAT( IX,  7E11.1*) 

2 

where  P3  was  just  the  constant  hjt/k  t)  .  If  KK  4  3>  either  PI  or  P2  is 
zero.  For  some  problems  with  aperture  excitation  very  near  the  poles  of  the 
body  of  revolution,  P2  was  negative.  This  apparent  error  is  due  not  to  the 
method  of  moments  itself,  but  to  the  approximations  made  in  obtaining  the 
source  and  field  integrals. 

Do  loop  9  calculates  complex  numbers  El  and  E2  whose  magnitude  is  the 
square  root  of  the  gain  by  pre-multiplying  the  matrix  of  current  coefficients 
TI  by  the  measurement  voltage  matrices  WR  previously  supplied  by  the  sub¬ 
routine  PLANE.  The  proper  normalization  is  obtained  in  statements  190  and 
191.  El  and  E2  at  1  and  NT  are  expected  to  be  zero.  Just  to  be  definite, 
the  b  statements  before  do  loop  139  assign  the  phase  0  to  El  and  E2  at  1 
and  NT.  The  magnitude  of  El  squared  is  the  gain  associated  with  the  t 
directed  aperture  electric  field  and  the  magnitude  of  E2  squared  is  the 
gain  associated  with  the  0  directed  aperture  electric  field.  Since  the 
measurement  voltage  submatrices  R^0  and  R^  are  zero,  the  t  directed  electric 
current  will  excite  only  a  9  directed  far  field  and  the  t p  directed  current 
only  a  p  directed  far  field.  Do  loop  139  converts  the  coefficients  TI  of  the 
current  expansion  functions  to  actual  currents  per  unit  length  by  dividing  by 
the  cylindrical  coordinate  radius  RH. 

After  do  loop  139>  the  remainder  of  the  program  is  devoted  to  printing 
and  plotting  the  electric  currents  and  gains.  The  currents  are  printed  in 
either  do  loop  l^o  or  do  loop  IU3.  These  currents  in  amperes  per  unit  length 
are  independent  of  <p  and  possess  real  and  imaginary  parts.  The  gains  are 
printed  in  either  do  loop  136,  139  or  l67>  In  the  heading,  GO  is  the  gain 
associated  with  the  t  directed  current  and  G 0  is  the  gain  associated  with  the 
0  directed  current.  E9  and  Bp  are  the  square  roots  of  the  respective  gains. 
ANG  E 0  and  ARC  Bp  are  indicative  of  the  phases  of  the  9  and  0  directed 
electric  fields  in  the  fur  zone. 

In  do  loop  171,  M  =>  1  obtains  plots  of  the  t  directed  current  JT  and 
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gain  GO  while  M  =  2  obtains  plots  of  the  <p  directed  current  J<p  and  gain  G0. 

Do  loop  172  finds  XI  and  X2,'  the  maximum  magnitudes  of  the  real  and  imaginary 
parts  of  the  current.  Do  loop  13  normalizes  XI  and  X2  to  full  scale.  K3 
and  k4  represent  the  real  and  imaginary  parts  of  the  currents  and  K1  and 
K2  their  corresponding  abscissas.  CALL  RE0RD(K1,K3,NM)  arranges  K3  in 
descending  order.  The  abscissas  K1  are  arranged  accordingly  to  maintain 
the  original  correspondence  between  K3  and  Kl.  CALL  RE0RD(K2,K4,NM)  operates 
similarly  upon  K2  and  K4.  The  current  is  plotted  in  do  loop  20.  One  line 
is  printed  for  each  J.  Data  is  plotted  on  a  given  line  only  if  the  ordinate 
of  the  data  is  greater  than  or  equal  to  the  ordinate  of  the  line.  If  some 
data  is  out  of  range,  that  which  is  too  large  accumulates  at  the  top  of  the 
plot  while  that  which  is  too  small  is  ignored.  Do  loop  80  prepares  the  gain 
to  be  plotted  on  a  logarithm  scale  via  K3  and  Kl.  CALL  RE0RD(K1,K3,NT)  arranges 
K3  in  descending  order  while  still  keeping  track  of  Kl.  The  gain  is  plotted 
by  do  loop  87.  At  the  end  of  the  program  execution  is  transferred  back  to 
read  statement  50.  If  there  is  another  set  of  punched  card  data,  the  pro¬ 
gram  recycles.  If  not,  execution  stops. 

Let  NP  be  either  the  original  NP  of  the  data  if  the  generating  curve 
does  not  close  upon  itself  or  2  plus  the  original  NP  of  the  data  if  the 
generating  curve  closes  upon  itself.  When  NP  >  43  or  NT  >  73>  some  vari¬ 
ables  require  larger  dimensions.  These  variables  are  listed  along  with  their 
minimum  dimensions. 

COMPLEX  Y(4*NM*NM),  WR(2*NT*NM),  Tl(2*NM), 

E3(2*NM),  El(NP),  El(NT),  E2(NP),  E2(NT) 

COMMON  R(NP-l),  ZS(NP-l),  SV(NP-l),  CV(NP-l), 

T(4*NM),  TR(4*NM) 

DIMENSION  RH(NP),  ZH(NP},  DH(NP-l),  TJ(NM),  THR(NT) 

DIMENSION  Kl(NM),  Kl(NT),  K2(NM),  K3(NM),  K5(NT),  K4(NM) 

NP-3 

Here,  as  in  the  main  program,  NM  =  — jp  .  The  variables  E1,E2,  Kl,  and  K3 
serve  a  dual  purpose.  The  common  statement  in  the  subroutine  PLANE  should 
be  exactly  as  it  is  in  the  main  program. 


hi 


VII.  BACKSCATTERING- ,  OBLIQUE  INCIDENCE 

This  program  computes  the  backscattered  field  and  monostatic  radar 
cross  section  for  a  conducting  body  of  revolution  excited  by  a  plane  wave 
from  an  arbitrary  direction  of  incidence.  The  current  is  also  found  but  not 
printed  out.  The  input  consists  primarily  of  NP  points  defining  the  body 
contour  plus  the  [y]  matrices  for  all  modes  n  satisfying  (53).  The  output 
consists  primarily  of  the  backscattered  field  and  monostatic  radar  cross 
section  patterns.  The  cross  section  patterns  are  graphed  by  a  printer  plot 
routine.  The  program  listing  plus  sample  input-output  data  are  given  in 
Appendix  D. 

The  subroutines  PLANE  and  REORD  are  exactly  the  same  as  the  ones  in 
the  program  concerning  scattering  for  axially  incident  plane  waves. 

Punched  card  data  is  read  early  in  the  main  program  according  to 

50  READ(1,51>  END  =  52)  KK,  HP,  NT,  BK 

51  FORMAT(3I3,  El M) 

READ(1,53)(RH(I),  I  =  1,  NP) 

READ(1,53)(ZH(I),  I  =  1,  HP) 

53  F0RMAT(l0F8.t) 

The  punched  card  data  for  the  backscattering  program  is  the  same  as  that 

for  the  program  concerning  scattering  for  axially  incident  plane  waves 

except  that  NN  is  replaced  by  KK.  For  oblique  incidence,  the  R  's  of  equa- 

r  1]  ** 

tions  (77)  and  (81)  of  the  report  are  considered  for  only  n  =  0,1,2*  *  *KX-1. 

Up  to  and  including  do  loop  8,  the  programs  for  backscattering  and  for  axial 

incidence  are  essentially  the  same. 

Do  loop  k2  initializes  El,  the  measurement  for  a  9  directed  transmitter 
and  a  9  directed  receiver,  and  E2,  the  measurement  for  a  0  directed  trans¬ 
mitter  and  a  0  directed  receiver.  When  both  transmitter  and  receiver  are 
located  at  0  «  0,  there  are  no  80  or  08  polarization  components.  For  the  96 

and  dub  polarizations,  R  Y  V  =  R  Y  V  so  that  the  measurement  must  be 

**n  -n  -n  n  n  n 

multiplied  by  an  extra  factor  2  if  n  4  0,  Inside  do  loop  to,  the  subroutine 


PLANE  is  called  upon  to  provide  WR,  the  receiver  voltage  corresponding  to 
R  for  n  =  KK-1.  From  the  KK^  record  of  data  set  6,  statement  127  reads 
the  admittance  matrix  Yn  for  the  same  (M  -l)^11  mode.  Each  L  of  do  loop  4l 
denotes  a  different  polar  angle  0.  For  transmitters  in  the  direction  0  from 
an  origin  near  the  body  of  revolution,  do  loop  3  puts  the  currents  in  E 3  and 
E4.  The  transmitter  voltage  matrices  necessary  for  computation  of  the  cur¬ 
rents  have  been  abstracted  from  the  receiver  voltage  matrices  WE  by  merely 
changing  the  sign  of  the  elements  corresponding  to  and  R^\  Do  loop  4-3 

adds  the  measurements  R  Y  V  to  El  and  E2.  Toward  the  end  of  do  loop  4l. 

n  n  n  c  ’ 

SIG  00  and  SIG  00  along  with  their  square  roots  MAG  S00  and  MAG  S00  are 

-  2 
printed.  SIG  00  and  SIG  00  are  echo  areas  per  square  wavelength  (o/k)  obtained 

by  squaring  the  product  of  El  or  E2  with  PI.  PI  is  the  constant  k2T]/27t^2 
defined  earlier  in  the  program  just  after  do  loop  1. 

The  fixed  point  ordinates  K3  and  K4  and  abscissas  K1  and  K2  appearing 
in  do  loop  80  facilitate  plotting  the  echo  areas  per  square  wavelength  on  a 
logarithm  scale.  Do  loop  8o  also  repeats  the  final  printing  of  the  echo 
areas  per  square  wavelength.  Statements  14  and  15  arrange  K3  and  K4  in 
descending  order  and  do  loop  87  plots  the  echo  areas  per  square  wavelength 
in  the  same  way  that  the  bistatic  scattering  cross  sections  were  plotted 
in  the  program  for  axially  incident  plane  waves. 

If  NP  >  43  or  NT  >  73>  most  variables  require  larger  dimensions.  The 
affected  variables  are  the  same  as  those  of  the  same  name  appearing  in  the 
program  for  axially  incident  plane  waves  except  E3  and  E4  which  must  be 
dimensioned  at  least  as  large  as  TI  was  in  the  program  for  axially  incident 
plane  waves. 
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APPENDIX  A.  PROGRAM  AND,  SAMPIE  INPUT-OUTPUT  DATA  FOR  SECTION  IV. 

//  (0034, EE, 6, 2) , »MAUTZ,40E* ,MSGLEVEl=l 

//  EXEC  FORTGCLG,PARM.FORT='MAP' 

/ /FORT . SYS  I N  DO  * 

SUBROUTINE  LINE0(LL,C) 

COMPLEX  C(1),ST0R,STU,ST,S 
DIMENSION  LR(58) 

DO  20  1=1, LL 
LR( I)=I 

20  CONTINUE 
M1  =  0 

DO  18  M=  1  ,LL 
K=M 

DO  2  I=M,LL 

K1=M1+I 

K2=M1+K 

IF<CABS<C<K1 ) )— CABSIC (K2 ) ) )  2,2,4 
4  K=  I 

2  CONTINUE 
LS=LR(M) 

LR ( M ) =LR ( K ) 

LR(K)=LS 

K2=M1+K 

STDR=C(K2) 

J1  =  0 

DO  7  J=1,LL 
K1=J1+K 
K2=J 1+M 
STO=C { K 1 ) 

C(K1)=C(K2) 

C I K2 ) =STO/STOR 
J1=J1+LL 
7  CONTINUE 
K 1 =M 1+M 
C ( K 1 ) = l , /STOR 
DO  11  I =1 ,LL 
IF(I-M)  12,11,12 
12  K1=M1+I 
ST=C(K1 ) 

C(K1 )=0. 

J1  =  0 

DO  10  J=1,LL 
K 1 = J 1+ 1 
K2=J1+M 

C(Kl)=C<Kl)-CU2)*ST 

J1=J1*U 

10  CONTINUE 

11  CONTINUE 
Hi=Ml*LL 

18  CONTINUE 
J1  =  0 

DO  9  J*  1  *Ll 
I F ( J-LR ( J ) )  14,8,14 
14  IRJ»LR(J) 

J2=(LRJ-1)*IL 

21  00  13  I«1,LL 
K2* J2*  I 
K1*JUI 
S*C(K2) 

CIK2)=C<K1 1 
CU1)  =  S 


45 


13  CONTINUE 

LR( J)=LR!LRJ> 

IR!LRJ)=LRJ 

IF ( J-LR ( J )  )  14,8,14 

8  J1=J1+LL 

9  CONTINUE 
RETURN 
E.NO 

COMPLEX  A3,A4,Z ( 1600 ) ,GS ( 40 ) , 0 ( 5292 ) ,U 
DIMENSION  RH<43) ,ZH<43) ,DH(43) ,TJ(20) 

DIMENSION  SV(42)*CV(42),ZS(42)»R(42)» ANG ( 40 ) , AC ( 40 ) »CSM ( 1 20 ) 
DIMENSION  TP<80),T<80),TR<fi0),JK(4) 

REWIND  6 
U=(0., 1. ) 

50  REA()(1,51,END=52)  NN,NP,NPHI  ,BK 

51  FORMAT ( 3 1  3, E 14. 7 ) 

READ! 1,53) (RH( I ) ,I=1,NP) 

READ! 1,53) (ZH< I ),I=1,NP) 

53  FORMAT ( 10F8 ,4 ) 

76  WR I TE ( 3, 54 )  NN,NP,NPHI , BK 

54  FORMAT ( IX // *  NN=',I3,'  NP=',I3,'  NPHI=',13,‘  BK=',E14.7) 

55  FORMAT (IX/1  RH* ) 

WRITE! 3,55) 

WRITE! 3,46) (RH(I), 1=1, NP) 

46  FORMAT! IX, 10F8.4! 

WRITE! 3,56) 

56  FORMAT! IX/'  ZH* ) 

WRITF(3,46)!ZH( I),I=1,NP) 

IF((RH(1)-RH(NP)).N£»0.«0R.(ZH(1)-ZH(NP))*NE»0.)  GO  TO  5» 

RH ( NP+ 1 )  =RH ( 2 ) 

ZH ( NP+1 )  =  ZH( 2  ) 

RH <  NP+  2 ) =RH<  3 ) 

ZH ( NP+2 ) =ZH ( 3 ) 

NP=NP+2 

5R  DO  57  U2,NP 
12=1-1 

RR1=RH( I )-RH( 12) 

RR2°ZH( I )-ZH! 12) 

OH!  12)=SORT(RRl*RRl4RR2*RR2) 

ZS ( I  2 ) = . 5*  t  ?H( I )4ZH!I2)  ) 

R! I  2 )  =  ,5* ( RH ( I )4RH! 12)  ) 

SV (  I  2 ) =RR 1 /OH { l 2 ) 

CV ! ! 2 ) =RR  2/OH (12) 

57  CONTINUE 
KG=NP-1 
NsK0/2 
NMaN-l 
NM2=NM*2 
NM4=NM*4 
NZ=NM2»NM2 
NGsKG»KG 
M5=NN>2 

M6=NN44 

FHaNN 

FM2=NN*NN 

P i °3. 141593 

FT  A  =  376,707 

DPaPl/NPHt 

CA=ex*8R*ETA 
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CQ=ET  A 
SS=0. 

DO  117  1  =  1, NM 
Il  =  2*(  I-l)+l 
12=11+1 

ss=ss+oh<  i  n+CHi  i2) 

TJ (  I  )  =  SS 
117  CONTINUE 

WRITE! 3,118) 

HR  FORMAT  <  1 X  /  *  TJ') 

WR I TE ( 3,46 )  (TJ!I),I=1,NM) 

DO  2  J=1 ,NPH I 
ANGl J ) = 1 J~,5 )*DP 
AC ( J ) =C0S ( ANG! J )  ) 

2  CONTINUE 
M3=0 

DO  10  MM=M5,M6 
M1=MM— 3 
M2=M3*NPHI 
00  11  K=1 ,NPHl 
K1=M2+K 

CSM ( K1 ) sOP^COS (M1*ANG ( K )  ) 

11  CONTINUE 
M3=M3+1 
10  CONTINUE 

DO  16  J=1,KG 
DEL=.5*DH(J) 

DEL 1=DH l J ) #BK 
AA=DP*R( J)*DEL*RK 
on  17  1=1, KG 
Z3=ZS( J)-ZS( I  ) 

RR1=SV! J)*RI J)+CV( J)*Z3 
RR2=-SV( J)*R(  I ) 

RR3=R 1 J )*R l J ) +R I  I ) *R { I )+Z3*Z3 
RR4=-2.*R( J)*R(  I ) 

X1=ABS!R(  J5-RU  ))+ABS!Z3) 

DO  5  K=1,NPH1 

IFlK.NE.l.OR.Xl.NE.O.)  GO  TO  7 
X=R ( J ) *0P 

XX=SQRTIDEL*OEL+X*X) 

W1=(X*AL0G! (OEL+XX ) /X )+D6L*AL0G( I X+XX ) /DEL ) ) /AA 
w2=-l . 

GS  (  l)=Kl+U*W2 
GO  TO  5 

7  Y=ABS«RRl+ACfK)*RR2) 

RD=RR3+RR4«ACIK) 

RK=8K*S0RT IRO) 
t)2  =  RD-Y*Y 
Y ) = Y-DEL 
Y2  =  Y-*0EL 

»1®S0RT( Yi*Yl+02) 

R2aSQRT(Y2*Y2+02) 

IF(Yl)  72,73,73 

72  T  1  N=  ALOr,  M  -Y 1 +R  l)  •  <  Y  2  ♦«?  )  /!)?  ) 

GO  TO  ?S 

73  tIN=ALOG(  (Y2+«?)/f  Y1*RU  > 

25  SN=S1NIRK) 

C$=COSUK) 

GS  t  K  )  =  ( C S~U*SN )  •  {  TlN-U»(HK»OM(  J)-RK*TIN)  )/t)EU 
5  CONTINUE 
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M3=< J-l )*KG+I 
00  A 8  MM* 1 » 3 
M1=MM-1 
M4=M  1  <«NPH  I 
M2=M1*NG+M3 
G ( M2 ) =0 . 

00  13  K=1,NPHI 
K?-=K+M4 

G( M2 1 =G( M2  I +GS  <  K ) *CSM  <K2 I 
13  CONTINUE: 

Afl  CONTINUE 
17  CONTINUE 
1A  CONTINUE 

00  74  J=1,NM 
J2=2*( J-1I+1 
J3= J2+1 
J4= J3+ 1 
J5=J4+1 
JA=4*( J-1I+1 
J7=jA+ 1 
JB= J7+ 1 
J9=J8+1 

OEL1=OH( J2I+DHI J3) 

0EL2=0H( J4J+0HI J5) 

TP( JA)=OH< J2I/0EL1 
TP(  J7)=0H< J3I/DEL1 
TP ( J8)=-UH{ J4I/0EL2 
TP( J9)=-0H( J5I/DEL2 
T( JA)=OH( J2!*0H( J2I/2./DEL1 
T(J7)=0H( J3)*(DH( J2I+0HI J3I/2.I/OE11 
T( Jfl)=0H(J4}*<0H| J5I+0HI J4) /2. I/DEL2 
T( J9)=0Ht J5I*DH(J5)/2./D?L? 

74  CONTINUE 

00  75  J=1,NM4 
TR( J)*T< J> 

75  CONTINUE 

115  !F(  (ZHU)-2H(NP-2)  ).EO.O..AND.(RH(1>-RHINP-2)).FO.O,)  GO  TO  7« 
1 F ( RH ( 1 ) )  77,23,77 
77  Otll=t)H(  1  J+OH(  ?) 

TM 1 )=QH( 1 )•( 1 .+(0H<2>+0H( 1 1/2. I /DEL  1) 

TR(  2 ) ■  f)H (  2 ) * t  l,+l)M(2l/2»/OELl  > 

23  IFIRHINPI)  79,7H,79 
79  J 1  =  ( N-2 I *4+3 
J2=Jl+l 

06L2»DMINP-2)+0M(KC,I 
TP  ( J1  )»DH!NP-2 1* 1 1  «  +  ()M(NP-2 1 /  2./OEL21 
11 A  TR I J2 1  «(>H(  KG  I  •  I  l.  +  (0H(NP-2)+0MUG)/2.  J/OE12  5 
7«  00  30  Ja 1 ,NM 
JL*< J~II*NM2 
J3* ( J-l ) *4 

ji=2*(j-n 

00  31  I = 1 ,NH 
L 1  a JL+ I 
l 2*L l+NM 
t3«NM*NM2+Ll 
L4=L3+NH 
2111 )*0. 

21121*0. 

2113) *0. 
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Z ( L4 ) =0, 

1 1  =  2*( 1-1) 

13=1 I-l)*4 
no  70  JJ=1,4 
J2=J1+JJ 
J7=J3+JJ 
00  71  11=1,4 
12=11+1! 

17=13+11 

J4=(J2-1)»KG+I2 

J5=J4+NG 

J6=J5+NG 

SS=SV (12) *SV !  J2 ) 

CC=CV( 12)*CV(J2) 

A3=.5*(G( J6)+G(J4) ) 

A4=.5*!G( JA1-GI J4) ) 

?  (L1)=Z(L1)  +  1CA’>T  ( 17)*T(J7)*C  SS*A3+CC*G(  J5)  )-Cvi*TP<  I7)*TP(  J7)*G(  J5 
1 )  )  »U 

Z!12)=ZCL2)+CA*SV(  J2)*TR(!7)*T!  J7)*A4-FM*CQ*G(J51*TRU7)*TP(  J71/P1 
112) 

Z(13)=Z(L3)-CA*SV(  12)=>T<  I  7  )  =»TR  <  J7 > *A4  +  FM*C0*G  (  J5 )*TP( ! 7 )*TR <  J7)/R( 
1J2) 

Z  (L4)=Z(L4)  +  {CA*A3-FM2«'C0/R(I2)/R(  J2)*G(  J5)  ) *TR  (1 7  )*TR  ( J7  ) *1) 

71  CONTINUE 
70  CONTINUE 
31  CONTINUE 
30  CONTINUE 
HI  CAU  LIN60(NM2,Z) 

WRITEtM(Z(n,l  =  l,NZ) 

JK  (  1  )  =  1 
JK(2)=N 

JK( 3)=NM2*N  (+1 
JM4)=JK(3)+NM 
1)0  93  J=l,4 
K  l  = JK ( J ) 

WR I TF ( 3, 24 )  J 
24  FORMAT (IX/1  Y  *  ,  1 1 ) 

00  92  1=1, NM 
K?nKl+NM-l 

<14  WRlTE(3»flH)(ZIK)*K=Kl*K2) 

HR  FORMAT!  1X.10GU.41 
K 1 =K 1 *NM2 
92  CONTINUE 
43  CONTINUE 
GO  TO  50 
52  STOP 
END 

/* 

//G0.FT04F001  00  0SNAHE=EEO034.REVl ,0!SP«OLO»UN! T=2314,  X 

U  V0LUME=SeR=SU0004,0CR=!«cCFM«v,BlKStZE=lft00.LRECL=1794) 

//gu.sysin  00  • 

001041  20  0. 4 AS 9995 E +00 


0.0 

0.0R6R 

0.1734 

0.2605 

0.3473 

0.4341 

0.5209 

0.6078 

0.6944 

0.7814 

0.8682 

0.9551 

1.0419 

1.12R7 

1.2155 

1.3024 

1 . 3892 

1.4740 

1.5428 

1.649? 

1.7365 

1  .0233 

1 .9101 

1.9970 

2.0R3R 

2.1706 

2.2574 

2.3442 

2.431 l 

2.5179 

2.6047 

2.6837 

2.6863 

2.5969 

2.4184 

2.1570 

1.R216 

1.42J8 

0.9772 

0.4971 

-0,0000 

0.0 

0.4924 

0.9R4R 

1.4772 

1.9696 

2.4620 

2.9544 

3.446^8 

3.0392 

4.4316 

4 .9240 

5.4144 

5.90RR 

6.4013 

6. 8937 

7 ,3861 

7.H7R5 

8.3709 

8.8633 

9.355? 

9.R481 

10.3405 

10.R329 

11.3253 

11*6177 

12.3101 

12.8025 

13.2949 

13.7873 

14.2797 

14.7721 

15.2657 

15.7650 

16.2562 

16.7225 

17.1478 

17.5177 

1 7. R 195 

18.0427 

18.1798 

IR.226Q 
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APPENDIX  B.  PROGRAM  AND  SAMPLE  INPUT-OUTPUT  DATA  FOR  SECTION  V. 

//  (0034, EE, A, 2),  *MAUTZ,J0E',MSGLEVEL=1 

//  EXEC  FOR TGCLG,P ARM. FORT  = 'MAP • 

//FORT . SYS  IN  OD  * 

SUBROUTINE  PL ANE ( VVR ,THR ,NT > 

COMPLEX  VVR( 1 ) ,A5,A6,U 

COMMON  U,R(42),ZS(42),SV(42) ,C V < 42 > ,BK ,NP ,NN, T ( RO ) ,TR(«0) 
DIMENSION  BJ( 126) ,THR < U ,FK ( 20) 

KG=NP-1 
NM=KG/2— 1 
M2=NN+2 

A5=2.*3.141593*U#*(NN+1) 

NV=NM*4 
FK( 11  =  1. 

DO  153  J=1,M2 
J1=J+1 

FK  { J 1 )  =FK  ( «J )  *  J 
153  CONTINUE 

DO  156  L=1 ,NT 
Ll=( L-l )*NV 
CS=COS( THR (L ) ) 

SN=SIN(THR(L)I 
BCS=BK*CS 
00  302  J=1,KG 
X=R( J)*BK*SN  . 

J1  =  J 
1 1=NN 

I F ( II)  303,304,303 

304  11=11+1 
Jl= Jl+KG 

303  DO  305  JJ=I1,M2 
I F ( X-l • E-5  )  1,1,2 

1  IF(JJ-l)  3,3,4 

3  B.l  (  J 1 )  =  1  • 

GO  TO  306 

4  BJ  (  J 1 1=0. 

GO  TO  306 

2  RH=X/2. 

RH2=RH*RH 
RH3=RM**( JJ-1 ) 

B J ( J1)=RH3/FK( JJ1 
SS=BJ ( J 1 ) 
ft  SST=SS*l.E-7 
DO  155  K= 1 , 20 
SS  — SS*RH2/K/(K  +  JJ-1) 

HJ  ( J 1 ) =RJ I J 1 ) +SS 
IF(ABS(SS)-SST)  306,306,155 
155  CONTINUE 
STOP  155 

306  J  1  ■ Jl+KG 

305  CONTINUE 
302  CONTINUE 

IF(NN)  307 ,308,307 
30fl  00  309  J= 1 »KG 
Jl=J+2*KG 
BJ(JI»-RJU1I 
309  CONTINUE 

307  00  300  J«l,NM 
J1=J+L 1 
J2=J1+NM 


J3=J2+NM 
J4= J3+NM 
VVR !  J 1 )  =0  • 

VVR(J2)=0. 

VVR ( J3 ) =0  . 

VVR(J4)=0. 

00  301  1=1,4 
I 1  =  2*( J-l )+I 
14=4*  <  j-u+i 

l 2= I 1+K6 
13= ! 2+KG 

AA=(COSi  ZS( II)* 8C S ) +U* S 1 N ( ZS( II) *BC S ) ) *A5 
BJ1= ( BJ( 1 3 )+BJ ( 1 1 ) )*.5 
BJ2=(BJ( 1 3 )-BJ ( 1 1 ) ) *. 5 

VVR(  J1)=VVR(  J1 ) *Aft* ( CS*SV  < Il)*BJ2+SN*CVn 1 ) *BJ (I  2 > *U ) #T ( 1 4 ) 
VVR(J2)=VVR(J2)+AA#CS*BJ1*U*TR(I4) 

VVR( J3)=VVR( J3)-AA*SV< II )*BJ1*U*T(I4> 

VVR  ( *14  )=VVR  I J4 )  +A6*BJ2*TR  ( 14 ) 

301  CONTINUE 
300  CONTINUE 
15 A  CONTINUE 
RETURN 
END 

SUBROUTINE  REORD! K 1 ,K3, L ) 

DIMENSION  K 1 ( 1 ) »K3 ( 1 ) 

DO  81  J=  1 ,  L 
K.fl=K3(J) 

KA=J 

DO  82  I = J  »L 
IF ( K3( I1-K8)  82,82,84 
84  K8=K3( I ) 

K  A=  I 

82  CONTINUE 

K3(KA)=K3( J  ) 

K3 ( J)=K8 
K8=K1 (KA) 

K1 IKA ) =K1 1  J ) 

Kl( J)=K8 
81  CONTINUE 
K3tl*l)=-1 
RETURN 
END 

COMPLEX  A3,Yt 1A00) ,VVR <58401 ,T! (40) , E 3 1  20) ♦ E 1 ( 73 ) , E2 ( 73 ) ,0 
COMMON  U,R(42),ZS(42>,SV(42),CVi42),BK,NP,NN,T(80),TR(P0) 
DIMENSION  RM(43),ZH(43),0HI42J,TJ(20) ,INT( 11),THR(73) 
DIMENSION  AA| 105),K1I73),K2(73),K3( T3) ,K4(73) 

DATA  AA(1) ,AA< 104) ,AAl 105)/*  ,,»X*,"G*/ 

DO  107  1=1,102 
10?  AA(  U1)=AA{  I  I 
U=(0.,1.  ) 

ETA=3?A.707 

ETA2=ETA*2. 

P!=3. 141593 
PR=180,/P! 

REWIND  A 

50  READ! 1 ,51,EN0=52)  NN,NP,NT,BK 

51  FORMAT!  3 1  3, El*,. 7 ) 

READ!  1 ,53HRMt  l),l  =  l.NP) 

REA0I1,S3MZH(  I),I=1,NP) 

S3  FORMAT ( 10FR.4) 


57 


76  WR I TE ( 3 , 54 )  NN,NP,NT,BK 

54  FORMATdX//'  NN=',I3,'  NP=',I3,'  NT='»I3,'  8K=»,E14.7) 

WR I TE ( 3 , 55 ) 

55  FORMAT (IX/1  RH'  ) 

WRITE(3,46) 1RH1 I),1=1»NP) 

46  FORMAT (lXtlOFft.4) 

WR ITE 1 3,56) 

56  FORMAT (IX/'  ZH» ) 

WRITE13.46) 1ZH(  I)  ,I=1,NP) 

KL=  1 

126  IF((RH(1)-RH(NP)).NE»O..OR.(ZH(1)-ZH(NP)).NE»0»)  r,0  TO  5« 
KL  =  0 

RH1NP+1  )=RH(  2) 

ZH1NP+1 ) =ZH( 2 ) 

RH ( NP+2 ) =RH 1 3 ) 

ZH { NP+2 ) =ZH ( 3 ) 

NP=NP+2 

5fl  DO  57  1=2, NP 
12=1-1 

RR1=RH( I )-RH ( 12) 

RR2=ZH( 1)-ZH( 12) 

1)H(  I2)=S0RT(RR1*RR1+RR2*RR2) 

ZS<  12)  =  .5*(ZHU  )  +Z H ( 1 2  ) ) 

R (  !2)=.5*(RH( I )+RH( 12) ) 

SV( I2)=RRl/OH{ 12) 

CV( 12)=RR2/OH{ 12) 

57  CONTINUE 

OT  =  P I / ( NT-1  ) 

DO  1  J=l,NT 
THRU)  =  I)T*  1  J-l  ) 

1  CONTINUE 
NM=(NP-3)/2 
NM4=NM*4 
NM2=NM*2 
NZ=NM2*NM2 
UO  74  J= 1 , NM 
J2=2* ( J-l ) +1 
J3=J2+1 
J4=J3+l 
J5= J4  +  1 
J6  =  4* ( J- 1 ) ♦  1 
j7=J6+l 
JR=J7+1 
J9=J8+l 

DEL  1 =UH ( J2 ) *DH (J3 ) 
l)EL2=OH(J4 )  *DH(  J5) 

T(  J6)=DHU2)*OH(J2  I/2./DEL1 
T(J7I=I)H|  J3)«(OHU2)+t)H(  J3)/2.  )/06ll 
T  ( J8)«0MU4)*10H(  J5)*OHU4»/2.)/OEt2 
TU9)=t)HU5)*l)H(J5)/2./D6L2 

74  CONTINUE 

00  75  J=1,NH4 
TR(J)oT(J> 

75  Ct)NTlNUE 

115  IF(KL.EO.Q)  GO  TO  78 
IFCHHim  77,23.77 

77  OEU  =  W l  l)«t)M(2) 

T  a |  1  )=OH( 1 )•{ l.«IDM(2t*0H| 1 )/2. )/DEL 1 ) 

TR  <2 )  =0M( 2) *(  1  .♦Oh (2 )  /2./DEU  ) 


23  IP ( RH( NP ) )  79,78,79 
79  Jl= ( NM-1 ) *4+3 
J2= Jl+1 

DEL2=DH ( NP-2 1 +DH ( NP-1 ) 

TR(Jl)=DH(NP-2)#( i ,+DH( NP-2 ) /2./DEL2 ) 

TRI J2)=DH< NP-1)* (1.+IDH1NP-2 l+DH (NP-1) /2.J/DEL2) 

78  SS=0. 

DO  7  1=1 ,NM 
11=2*1 1-11+1 
12=11+1 

SS=SS+UH( 1 1 )+OH( 12) 

T J ( I >  =SS 

7  CONTINUE 
OEL=TJ(NM) 

IF(KL.NE.O)  DEL=OEL+OH ( NP-2 ) +DH ( NP— 1 ) 

OEL=OEL/10. 

DO  8  J  =  1  *  NM 
TJ(J)=TJ( Jl/DEL 

8  CONTINUE 

85  CALL  PLANE(VVR,THR,NT) 

127  REA0(6MY(  I  ) , 1= 1 ,NZ  ) 

DO  108  INC=1  ,2 
J3  =  0 

1PI INC.EO.l)  J3=NM4*(NT-1> 

DO  2  J=1,NM2 
T I ( J ) =0 . 

DO  3  1=1, NM 
J1=J+(I-1)*NM2 
J2=J1+NM*NM2 
l 1= I ♦ J3 
I  2= 1 1+NM 

T I ( J ) =T I { J ) -Y ( J1 ) *VVR ( 11 ) +Y  C  J2 1 *VVR  ( I  2 1 
3  CONTINUE 
2  CONTINUE 
DO  9  J  =  1 ,  NT 
SllJ)=0, 

E2<  J ) =0. 

Jl=( J-1)*NM4 
DO  10  1=1, NM2 
11=01+1 
12=1 1+NM2 

eiiJi=ei(j)+vvRini*Tim 

E2U)=P2(  J)+VVRU2)*TK1 ) 

10  CONTINUE 

9  CONTINUE 

J5=t  2-lNC)»(NT-l 1  +  1 

A3  =  CARS(E1U5)  1/E11  J51*(RK»BK*ETA/2,/P|/SORT(Pl  1  1 
00  11  J=1,NT 
EK  J»=Elt  J1+A3 
E21 J 1  =E21 .1)  *U*A3 

11  CONTINUE 
WUITE13.1 101 

110  E0«MAT('t*,2X,'T«,<.X,*REAL  JT*,lX.*fMAG  JT', 2X,*MA.-  JT«,lX,»«fAL  J 
10',1X,' 1MAG  JO*,?X«,MAG  JO') 

MR  1 Tp ( 3 , 109 1 

109  FORMAT t  »♦* ,37X,'/« ,TX, •/• ,7X,'/'  ) 

DO  128  Ja l ,NM 

J 1  a J*NM 

J2  =  2«M  J-ll+3 

TIUlaTH  J)*ETA2/UMU2) 


E3I J)*TI < J1)*U#ETA2/RH(J2) 

12ft  CONTINUE 

011  129  J=l, NH 
J 1 = J+NM 
J3= J-l 
J5=J+1 

I F  t J.NE. l.AND.J.NE.NM)  GO  TO  125 

J3=J 

J5  =  J 

IF(KL.EQ.l)  GO  TO  125 

J3  =  NM 

J5=J+1 

IFiJ.fcO.il  GO  TO  125 

J3=J-1 

J5=l 

125  TI  (J1  )=.25*<E3I J3 1 +2 .#E3 ( J 1 +E3 ( J5) 1 
129  CONTINUE 

00  4  J=1 ,NM 
J 1  =  J  +  NM 
X2=RE  AL  ITKJ)) 

X3=A I  MAG ( T I  I J  1  ) 

X4=CABS  ITHJI) 

X5=REAL<  TI  (  J 1 )  ) 

X6=AIMAG(TI< Jl) ) 

X7=CABS I T I ( J1 )  ) 

WRITE (3. 124)  TJI J) ,X2,X3,X4,X5,XA,X7 
124  FORMATIlX.fc5.2t6F 8. 3) 

4  CONTINUE 
WRI Ib(3.112) 

112  FORMAT  I  •  1 1  »  '  0',4X,*SIf>  X0».2X,»S!G  X0',2X,»MAG  SXO'tlX.'ANG  9X0 

1, IX, 'MAG  SXO'.lX.'ANG  SXO' ,  IX,  M.SIG  X0‘  ,  IX,  ‘UK,  XO') 

WRITE(3,113> 

113  FORMAT! •+• ,2X, »-',9X, •- * , 7X , ' / • , «X , •- » , 7X , •- » , 7X, • / • , 7X , • / • , ?X , • - 
1.7X, »/« ) 

00  12  J=1,NT 
X1=THR( J)*PR 
X4=CAHS( El( J  > ) 

Xft=CAHS(fc2t Jl) 

X2=X4*X4 

X3=X6*Xft 

X5=PR*ATAN?<  AIMAt,(El<  J)  l.RFAUFl!  Jl!) 

X7=PR*ATAN2t AIMAG«E2( Jl) ,R6Al!E2!  Jl)  1 
Xft= AL0G1 0(  X2 ) 

X9  =  ALOG 1 0  (  X 3 ) 

WRITE (3, 111)  XI ,X2,X3,X4,X5,X6,X7,XH,X9 
111  FORMAT! lX,F5.l,3Ffl.3,F8.1,F«.3,FR.l,2FR.31 

12  CONTINUE 

00  13  J=1,NH 
KH  J)bTJ(  J)*10.+3.5 
K2( J)=K1 (J) 

K3!J)-CARS(TMJH*10.>.5 
Jl sJ  +  NM 

K4( J)=CABS(TI( Jl ) 1*10. *.5 

13  CONTINUE 

14  CALL  REO«0(K1,K3,NM) 

15  CALL  RFnR()U?,K4.NM) 

00  104  J=l,ll 

INTI J)=J-1 
104  CONTINUE 
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K  =  5 
K5  =  1 
K6=l 

WR  I  TE ( 3, 106 ) 

10ft  FORMAT ( • 1» ) 

00  20  J=l»51 
J 1=5 1-J 
WR I TE ( 3 , 25  ) 

25  FORMAT ( '  I  •  ,<m,'  I ' ) 

IF! < J-l)/5*5-( J-l >)  21 ,22,21 
22  WRITE(3,123) 

123  FORMAT(  •+•  ,3X,  •  —  •  ,97X,  •  —  •  ) 

! F (  (  J-l)/10*10~(  J-IM  21,122,21 
122  WR I TE ( 3, 24 )  K 
24  FORMAT( '+• ,12) 

K=K-1 

IF(J.NE.l)  GO  TO  21 
WR  I  TE  t  3,  1 16  ) 

lift  FORMAT! '+• ,4X,50( ' — • ) ) 

WRITE! 3,47) 

47  FORMAT! •  +  » ,BX,19i • I  * ,4X) ) 

21  I F ( K3( K5 ) . LT  «  J 1 )  GO  TO  26 

60  K8=K 1 ( K5 ) 

WRITE (3,48) (AA(I) , I =1 ,K8 ) , AA( 104) 

48  FORMAT! '+*,10541) 

K5=K5+1 

I F ( K3 ( K5 ) ,GE , J 1 )  GO  TO  60 

26  IF(K4(K6),LT. Jl )  GO  TO  20 

61  K  H=K  2 ( K6 ) 

WRITE (3, 48) (AA( I ) , I = 1 ,K8 ) , AA < 105) 

K6=K  6+1 

I F ( K4 ( K6) .GE ♦ J 1 )  GO  TO  61 
20  CONTINUE 
WR I TE ( 3,47 ) 

WRITE(3,116) 

WR!T6(3,63)(INT!J),J=l,ll) 

63  FORMAT (3X,11( 12, 8X) ) 

WR I TF ( 3*  6A  ) 

64  FORMAT! 'O' ,20X,'X  X  X  PLOT  OF  MAGNITUDE  OF  T  DIRECTED  CURRENT  VER 
1SUS  LENGTH  T') 

WR 1 T  fc ( 3  ,65 ) 

65  FORMAT (21X, '0  0  O  PLOT  OF  MAGNITUDE  OF  0  DIRECTED  CURRENT  VFRSUS 
1L6NGTH  T • ) 

WRITE (3,66) 

66  FORMAT!  *♦•  ,48X,'/M 
DO  80  J= 1 , NT 

K1 ( J)=THR( J)*7?./PI+8.5 
K2( J)=Kl I J ) 

K3( J)=20.*AL0G10(CAMS(FH J) ) >+?0.S 
K4 ( J ) =  20. • ALOGIO (CAftS ( E2 ( J ) ) >+20.5 
80  CONTINUE 

16  CALL  REUROtKl ,K3,NT  ) 

17  CALL  RF0R0!K2,KA,NT  ) 

DO  105  J= l , 5 
InT(J)=! J-l)*45 

105  CONTINUE 
Xl=1000. 

*5=1 
K6  =  l 

WRITE! 3, 106) 
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DO  87  J=l,51 
J  1=51- J 
WR  I  T  b ( 3  *  88 ) 

88  FnRMAT!9X,*I',71X,«I') 

IF (( J-l )/ 10*10- ( J-l ) )92»90»92 

90  WRITE!3,9l)  XI 

91  FORMAT  (  •  +  •  »F7.?»  •  —  •  */S9X  *  •  —  •  ) 

X1=X1/10. 

IF(J.NE.l)  GO  TO  9? 

WRITE<3,93> 

93  FORMAT ( *  + ' ♦ 17X *7  I ' I ' tflX  )  ) 

WR I TF ( 3* 97  ) 

97  FORMAT! • ♦ ' ,8X , 73 ( « - *  )  ) 

92  IF!K3!K5).IT.J1>  GO  TO  94 
95  KB=K 1 ( K5 ) 

WRITE (3,48) ( AA { I ) , I=1,K8),AA( 104) 

K5=K5+1 

IF(K3(K5).GE.J1 )  GO  TO  95 

94  IF!K4(K4).LT.J1)  GO  TO  87 
9 4  K8=K2!K4) 

WRITE (3, 48)  (AAU  )  ,I  =  1,K8), AA!  105) 

K4=K4+1 

I  F I K4(K4 ) .GE . J 1 )  GO  TO  94 
87  continue 
WR I T  E ( 3  *93 ) 

WRITE l 3, 97) 

WRITE(3,98)  UNT!  J),J  =  1 ,5) 

98  FORMAT (7X, 3! I  3, 15X) . 14, 15X, 13) 

WRITEI3.99) 

99  FORMAT! *0' ♦ 15X, 'X  X  X  PLOT  OF  SIGMA  XO  OVER  LAM«OA  SOUAREO  VERSOS 
1  THETA* ) 

WKITEOtlOl) 

101  FORMAT!  '  +  ' ♦  3  7  X  ,  *  —  * ) 

WRITEI3.100) 

100  FORMAT! 14X,*0  0  0  PLOT  OF  SIGMA  XO  OVER  L AMROA  SOUAREO  VERSOS  THE 
IT  A 1  ) 

WRITEI3.102) 

102  FORMAT ( •  «  37X  , •/* ) 

WR 1 T  E ( 3 1 104 ) 

108  CONTINUE 
GO  TO  50 
52  STOP 
END 

/• 
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//  !0034,EE,6,2)  ,'MAUTZ,J()E'  ,MSGL6VEL  =  1 

//  EXEC  Fl)kTGCLG»P  ARM.  FORT*' MAP* 

//FORT . SYS  I N  DD  * 

SUBROUTINE  PLANE! VVR, THR, NT ) 

C (IMPLEX  VVR  (  1 )  *A5, A6,U 

COMMON  U,R(42I,ZS<42>  ,SV<42 ) ,CV(42  )  , BK ,NP, T ( RO > , TR ( 80 ) 
DIMENSION  BJI84) ,THR( 1 ) 

KG=NP-1 

NM=KG/2-l 

A5=2.*3. 1415<J3*U 

NV=NM*2 

DO  I  b A  L  =  1 »  NT 

L 1  =  < L-l )*NV 

CS=COS(THR(L) ) 

SN=S I N (  THR(L) ) 

BCS=BK*CS 
DO  302  J=1,KG 
J1=J 

X=R( J)*BK#SN 
UO  305  J J=1 «  2 
I F ( X-l .E-5 )  1,1,2 
1  IF(JJ-l)  3,3,4 

3  8  J  <  J  1 )  =  1  . 

GO  TO  306 

4  BJ(J1)=0. 

GO  TO  306 

?  RH=X/2 
RH2=RH*RH 
RH3  =  RH**M  JJ-1 ) 

8.1!  J1  )=RH3 
SS  =  8 J ( J l  ) 

H  SST  =  SS<>l.E-7 

00  155  Kb  1 , 20 
SS=-SS*RH2/K/(K+JJ-1 ) 

8 J ( J  1  )  =8,1 1  J  1  )*SS 
I F ( A8S ( SS I-SST )  306,306,155 

155  CONTINUE 
STOP  155 

306  J 1 = J 1 ♦KG 
305  CONTINUE 
302  CONTINUE 

DO  300  J=l,NM 
J 1  =  JH  1 
J  2*  J 1 ♦NM 
VVR! 01 ) =0. 

VVR ( J? ) =0, 

no  3oi  1=1,4 

I  1=2®! J-l )♦! 

J4«h«U-1  1*1 
l 2* I l*KG 

66* (COS I /St  I 1 1 *RCS I ♦U*SI N( Z S I 1 1» •RCS) )»A5 

VVRU1  )  *  VVR  I  J 1  )  ♦AA*  (C  S*SV(  I  1  )*HJ(  12USN*CV(  t  1  )*8J(  I  1  >*U)*T1  t  <• ) 
VVR(J2)=VVR(J2l*A6*8J( I2>*TR(I4J 
301  CONTINUE 
300  CONTINUE 

156  CONTINUE 
RETURN 
END 

SUBROUTINE  RpORD (K 1 ,K3 ,L  1 
01  HEN  SI  ON  Kim,K3UI 
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on  fli  j=i,i 

K8=K3< J) 

K4=J 

.  1)0  82  I  =  J,L 

I F ( K3 ( I ) -K8 )  82*82 *84 
84  Kfl=K3 ( I) 

K  4=  I 

82  CONTINUE 

K3 ( K4 ) =K3 ( J  ) 

K3( J)=K8 
K8=K 1 ( K6 ) 

K1(K4)=K1(J) 

K 1 ( J ) =K8 
81  CUNTINUE 
K3  ( L  +  l )  =-l 
RETURN 
END 

COMPLEX  A3,Y(1400),VVR(2920),TI ( 40 ) * E3 ( 40 ) , E 1 ( 7 3 ) , E2 ( 73 ) , U 
COMMON  U,  R( 42 ) »ZS ( 42 ) »  SV ( 42 ) »CV(42)*BK*NP,T(80) *TR ( 80 ) 
DIMENSION  RH(43),ZH(43),DH(42),TJ<20) » INT <11 ) , THR ( 73) 
DIMENSION  AA(llO) ,K1 (73) ,K2(20) ,K3 < 73 ) ,K4 ( 20) 

DATA  AA( l),AA(109),AA(110)/«  *,«X',*0*/ 

00  107  1=1,107 
107  AA ( 1+1 )=AA( I ) 

U=(0.,W> 

ETA=37A.707 
PI=3. 141593 
PR= 1 R0« / P 1 
REWIND  4 

50  READ! 1,51, EN0=52)  KK,NP,NT,BK 

51  FORMAT { 3  I  3 » E 14 ,7 ) 

READ!  l,53HRH(  I  ),I  =  1,NP) 

REAOt 1,53) (ZH( I ) ,1=1 *NP) 

53  FORMAT ( 10ER.4) 

WR I T  E ( 3 ,54 )  KK.NP.NT.BK 

54  FORMAT ( 1 X// 1  KK=',13,'  NP=',13,'  NT=',I3,'  BK=',E14.7) 
WR1TF(3,55) 

55  FORMATUX/*  RH*  ) 

WRITE(3,44)  (RMU)  ,IM.NP) 

44  FORMAT  I  IX, 10F8.4) 

WR  ITH3.54) 

54  FORMAT (IX/*  Z  H  •  ) 

WRlTE(3,44ICZHm,|al,NP) 
no  144  J=l,NP 
til J)=0. 

E2  ( J ) =0. 

144  CONTINUE 

IF(KK.FO.J)  GO  TO  40 
REAOt 1,148)161(1), 1  =  1 »NP  ) 

1 48  FORMAT (701). 4) 
w«ng(3.131) 

131  FORMAT! IX/'  60'  ) 

WR IT  t ( 3. 1 32 ) 

132  FORHATl -• ) 

WRITE  (3.  lARUtlH)  ,l*l,NP) 

149  FORMAT! IX, 7E1 1.4) 

IF(KK.EO.l)  GO  TO  41 

40  REALM  1,148)162(1), i«l, NR) 

WR|Tfc(3,13)> 
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WR ITE( 3* 133) 

133  FORMAT* •  +  /*) 

WRITE(3,169HE2(I),I  =  1,NP) 

41  KL=  l 

126  IFH  RH( 1 J-RHINP )).NE.O..OR.(ZH(l)-ZH(NP)).NE.O.)  GO  TO  5fl 
KL=0 

RH( NP+1 ) =  RH{ 2 ) 

ZH (NP  +  1 ) =ZH ( 2 ) 

El  (  NP+ 1 ) =E 1  (  2  ) 

E2 ( NP  +  1 ) =E2  (  ? ) 

RH( NP  +  2 ) =RH( 3 ) 

ZH(NP+2 )=ZH(3 ) 

£1 ( NP+2 )=E1  (  3) 

E2 ( NP+2 )=E2(3 I 
NP=NP+2 

58  00  57  I=2.NP 
12=1-1 

RR  1  =  RH (  D-RHU2) 

RR2  =  ZH{ I )-ZH(  1  2 ) 

OH ( I2)=SQRT(RR1*RR1+RR2*RR2) 

2  S  (  I2)=.5*(2H< D+ZHU2)) 

R( 1 2 ) = . 5* l RH( I ) +RH ( I  2  I  I 
SV{ !2)=RR1/0H( 12) 

CV( I2)=RR2/OH(12) 

57  CONTINUE 

OT  =  P 1 / ( NT- 1 ) 

00  1  J=1,NT 
THR  (  J  )  =DT t  J — 1 1 
1  CONTINUE 
NM=  t  NP-3) /2 
NM4  =  NM#>4 
NM?=NM*2 
NZ=NM2*NM2 
00  74  J=  1  »NM 
J2=2*(  J-l >♦1 
J3= J2+ 1 
J4= J3+  1 
J5= J4+ 1 
J6«4* ( J-l 1 ♦! 

J7=J6+1 
Jfl= J7+ 1 
J9=J«+l 

OEU  =  0M(  J? ) ♦OH t  J3 ) 

DEl2=0Htj4)>0H( J5> 

T*  J6)=0H«  J2>#OHU2)/2./OEU 
T  I J7 )  a  OH  (  J3 )  *  ( OH  ( J2 1  ♦OHf  J3)/2*l/t)El.l 
TUR)»OM(  J4>*(0HE J5 » ♦0M4  J4)/2,)/0El2 
T(J9l»0H(  J5)»t)H(  J5I/2./0E12 
7*.  CONTINUE 

00  75  J*1,NN4 
TR(J)aTU) 

73  CONTINUE 

115  IF! Kl, £0.0*1  GO  TO  78 
IF(RH( 11)77,23,77 
77  0EL1>0H( 11+0HI21 

t  r  1 1 1  *oh  i  ii*u.*ioH»2»-*OMcn/2.i/oeui 
TRt2)>0HI 21* ( l.*OH(2)/2./OEU  > 

23  IFIRH(NP) >79,78,79 
79  J1*(NN-1 1*4+3 
J?=J l+l 
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DEL2=DH(NP-2)+DH(NP-l ) 
TRUl)=OH(NP-2)>Ml.+OH<NP-2)/2./DEL2> 
TRU2)=DH<NP-l)*{l.  +  (OH<NP-2)+DH(NP-n/2.  )/DEL2) 
78  SS=0. 

DO  7  I =1 ,NM 
1 1  =  2* (  I-1I+1 
12=11+1 

SS=SS+DH (11) +OH (12) 

TJ ( 1 >=SS 

7  CONTINUE 
DEL=TJ(NM) 

IFIKL.NE.O.)  D£L=DEL+DH( NP-2 ) +DH (NP-1 ) 
t)EL=DEL/10. 

00  8  J=l,NM 
T  J ( J )  =  TJ ( J ) /DEL 

8  CONTINUE 

DO  44  J=1,NM 
J1=J+NM 
J2=2*  C  J-l ) 

J3=4*l J-l ) 

£3 ( J ) =0, 

E3( J1 1=0. 

00  45  1=1,4 
I 1=J2+I 
12=11+1 
I3=J3+I 

t3(J)=E3(J)-(ElUl)f61U2))*TU3) 

E3(J1)=E3(J1)-(E2(I1)+E2(I2))*TR(!3) 

45  CUNTINUE 

E3( J ) =63 ( J)*Pl 
E3IJ1 )=E3(J1)*PI 
4 4  CONTINUE 

185  RtAOU)(Ym,5  =  l,N7) 

P1=0. 

P2  =  0. 

P3=4.oPl/BR/BK/fcTA 
00  134  J=1,NM 
J l = J+NM 
T I t  J ) =0. 

T 1 1411=0. 

00  135  1=1, NM 
l 2= ( I-11+NP-2+J 
I3=12+<nm2M)*nm 
| 4o I *NH 

Til  J)  =  TI I  JI+YI I2J+E3I t ) 

TllJl)  =  TIun+YI!3»*E3(I4j 
135  CONTINUE 

P  1  =P  1  +C0N  Jf»  (  T  I  l.l>)*E3U> 

P2*P2+C0NJGin  (  Ji  I  I+E3I  J1 1 
134  CONTINUE 

WRITE  I  3, 184 j  Pl.P? 

184  EORMATI //1X.7E11 .4) 

P1  =  SCP.T|  A8SIPU+P3) 

P2  =  S2«TURSIP2»«P31 
IElPl.Nc.0.1  Pl+I./Pl 
IEIP2.NE.0.I  P2«i./P 2 
CALL  PLANE  I VVR, tHft,NT ) 

00  <»  .1=1, NT 
Eli J»=Q. 
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E2( J)=0. 

Jl= ( J-l ) #NM2 
00  10  1=1, NM 
I 1=J 1+ I 
12=1 1+NM 
13=1 +NM 

El ( J )=E1 ( J )+VVR! 1 1 )*TI ( I ) 

E2( J)=E2 ! J )+VVR< I21*TI ( 1 3) 

10  CONTINUE 

190  61 ( J ) =E 1! J  )*P1 

191  E2 ( J ) =62 ( J )  *P2 
9  CONTINUE 

fcllll=l.E-10 
E2( 1 ) = 1 . E-10 
El ( NT ) = 1 «E-10 
E2 ( NT )  =  1 « E-10 
DO  139  J=1,NM 
J?=2*( J-l 1+3 
J1=J+NM 

T1(J)=TI (J)/RH(J2) 

TI(J1)=TI(JI)/RH(J2) 

139  CONTINUE 

GO  T0( 145,146,137) ,KK 

145  J1=0 

WR I TE ! 3  »  13R  ) 

138  FORMAT! • 1' ,2X, *T',5X, 'REAL  JT»,4X,'!MAG  JT'J 
GO  TO  147 

146  J 1 =NM 
WR1TE(3,148) 

148  FORMAT  {  •  1  •  ,  2X,  1  T '  »  5X,  'REAL  JO'.AX.MMAG  JO') 

WR1TEI3.149) 

149  FORMAT! • ♦ • , 14X , • / • , 10X , • / • ) 

147  00  140  J=l,NM 
J2=J1+J 

WRITE  13.124)  TJ!J),TI(J2) 

124  F0RMAT(1X,F5. 2,4611. 3) 

140  CONT INUE 
GO  TO  150 

137  WRITfc(3,110) 

110  FORMAT (*1',2X,'T,,5X*'REAL  JT • , 4X , M MAG  JT',4X,'R6AL  JQ',4X,'!MAG 
1J0*  ) 

WR l TE ( 3, 109  J 

109  FORMAT! • ♦ * , 36X , • / • , 10X , « / * ) 

00  143  J= 1 »NM 
J  1  ® J  >NM 

WRIT613.124)  TJ!J).TI!J),TI(Jl) 

14J  CONTINUE 

150  GO  TO  ( 1 5 1 « 1 52* 1 53 ) ,KK 

151  WRITEI3.155) 

155  FORMAT ( *  1  0'.6X,'G0*,4X,'E0',4X,'ANG  60') 

WR I TE ! 3 » 1 54 ) 

154  FORMAT ! • ♦  -♦,7X,«~»,7X,'-*,9X,'-') 

00  156  J»l,NT 
X3»CAHS! E l( J ) ) 

Xl*THR(  J)*t»R 

X2sx3»x3 

X4*PR*AlAN2lA|HAC,ltl  i  J)  ), REAL  (El  !  J )  »  ) 

WRtT6!3,nn  Xl.X2.X3.X4 
1 t l  FORMAT! IX.F5.1,2F«.3,F8.1,2F8.3.FR.H 

156  CONTINUE 


- 
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GO  TO  157 
152  WR I TE ( 3t 155 ) 

•IR1TE(3*158) 

158  FORMAT ( *+  / * » 7X , * / 1 t 7X,  /  ♦ 9X»  I  ) 

00  157  J=1.NT 

X3=CABS(E2<J) > 

Xl=THR( J)*PR 
X2=X3<“X3 

X4=PR«ATAN2( A!MAG(E2( J  ) )  *RE AL t E2 l J) ) ) 

WR11H  3f  111)  X1*X2«X3»X4 

159  CONTINUE 
GO  TO  157 

1E0 '  ) 

WR  I  TF(  3f  161 )  |  .  |  -j  y  t/i  qv  #  i  / t  \ 

161  FORMAT ( 1 ♦  -•,7X,«-*.7X,*-»*9X.»-'.5Xt'/**7X,  /  ,9X, 

DO  167  J= 1 *NT 
X3=CAMS(E1( J) ) 

X  1  =  THR I J )  *PR 
Xj)  =  X3^X3 

X4=PR*ATAN2(  AIMAGI6K  J) )»  REAL  I  El  t  J> ) ) 

X6=CABS(E2(  on 

XTspR’J'AT  AN2  (AlMAGIE2(v))l  »REALI62(J)  )  ) 

WR  ITE(3« 111)  XI  ,X2t X3,X4»XS»X6,X7 
167  CONTINUE 
157  M 1 3 1 
N?  =  2 

IF(KK.EO.l)  m?=i 
IFIKK.E0.2)  Ml=2 
DO  171  M=M1,M2 
M3® ( M-l ) »NM 

X1  =  ARSIREAU  TI  (M3+1  )  )  ) 

X2®ABS(  AIMAC.m  (M3*l )  >  ) 

00  17?  J ® 1 » NM 

X3®AH5(R6Al(TI{Jt>n 
XAaABSI  A1MAGIT  I  I  JU  )» 

IF<  IX3-XD.GT.O.)  X1»X3 
IFUX6-X2).C»T,0.)  X2®X4 
172  CONTINUE 

DO  13  J*ltNM 
J1=J*H3 

KM  J)=T  jumo^A.5 

K2I  J)aKU  J)  _  * 

H3IJ)S25.6REALITI I  Jl ) )/Xl*?5»5 

K4U)s?5.«A|MAMTl(  Jl  I) / *2*25. 5 
13  CONTINUE 

CAU  «EUR0IK1,K3.NH) 

CALL  Rfcqf^lRZ.KA.NM) 

00  10*  J®  I » 1 1 
INTI 

10*  CONTINUE 
X  1  =  1  . 

*5=1 

*6=1 

WRITE  I  3* 10** 

106  FORMAT l ‘1*1 
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on  2C  4=1,51 

Jl=51-J  » 

WR  J  TE ( 3, 25 ) 

25  FORMAT ( 9X  »  • I • ,99X,  *  I  •  ) 

!H (J-l)/5*5-(J-l) )  21*22*21 
22  WR 1 TE ( 3 , 1 23  I 

123  FORMAT  (  •  ♦ '  ,  RX »  '  —  •  ,  97X,  •  —  •  ) 

122  WRITE! 3*29)  XI 

24  FORMAT  (•■*■•  ,F7.1  i  j 

Xl=Xl-.2 

IF  (  J.NE.  1 )  GO  TO  173 
HR  1 T  F ( 3,1141 

11A  FORMAT ( •  +  • ,9X,50( ' — • ) ) 

WR I T  E ( 3 ,47  I 

47  FORMAT (,+,,13X*19('I'*4X)) 

173  1  F  (  J.NE. 24)  GO  TO  21  j 

WR  1  TF  (  3 , 1 1 A ) 

21  I F ( K 3 ( Kb ) »LT  . J l )  00  TO  2A 
AO  Kfl  =  Kl ( Kb  ) 

WRITE! 3,4  ft) C  A  A ( 1 ) , 1  =  l , KR ) , AA ( 109 ) 

48  FORMAT ('♦•* 1 10A1 ) 

K5=Kb*l 

IF { K3! K5 ) «  GE • 41 )  00  TO  AO  j 

2A  IF(K4!KA).LT.J1>  GO  TO  20  1 

A1  KH=K2!  KA ) 

WRlTEI3»4fl)!AA!I),3=l»KR)*AA(110) 

K A=K A* 1 

IF(K4(KA).GE.J1)  00  TO  A1 
20  CONTINUE 

WRITE! 3*47 )  ’ 

WR  I  TE  I  3 •  1 1  A )  ! 

WRITEI3.A3H |NT( 41.4=1,111 
A3  FORMAT ( 8X , 1 1 ! 12,8X1/, IX) 

WRITE! 3, 174) 

174  FORMAT (32X,»X  X  X  PLOT  OF  (MEAL  J  I /MAX • 'REAL  J  ••  VERSUS  LENOT 
1M  T  •  I 

IF(M.EQ.l)  WR I TF I  3, 1 75  I 

175  FORMAT! ,54X, *1 • ,1?X, *T ■ )  | 

IHM.FO.?)  WR I TF ( 3, 1 7A ) 

17A  FORMAT!  ,54X  1 2X  ,  *0  54X  ,»/•,  1  2X  ,*/•  1 

UR  IT  F ( 3 , l HO ) 

IRO  FORMAT!  32*. •«  0  0  PLOT  OF  (IMAG  J  )/MAX‘'|HAf.  J  ••  VF«SWS  LF  NO  I 
|M  !'! 

IFIM.FQ.ll  WRI1F!3,175> 

IFCM.tQ.?)  WRITFI3.17A  | 

1A3  00  BO  J«l,NT 

Ml  4>=lMR!  J»*72./PI*».„ 

IFIM.FQ.il  K3IJ)»20.»A!.0f*10  SCABS  <F  t (J)  1 1*30.5 
|F(M.F0.2)  K3  i  4  1*20.  •AL'TO  1 0!  CARS  !  E2  <  J  J  >>*30.5 
BO  CONUNUF 

CALL  RF0RDU1  ,K3,NT  1 

00  105  4*1,5  | 

|NT!41=I4-I )**5 

105  Continue 
xi» ioo. 

K5*l 

Wfi  t  T F ! 3, IOA) 
nn  h7  4*1,51 

4 1 «  5 1  -  J  , 

WRITF13. R«>  ' 


! 
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«H  FlIRMATm,  '  l  * ,  7 1 X  »  '  !'  ) 

I F  (  (.1-1  )/ 10*10-1  J-l)  )92,<JO,<J2 
<10  WR  |  T  F  (  3 , 9 1  )  XI 

91  FORMAT ( ,F7.3* •  --'.ARX,'--*) 

X1=X1/10. 

lFU.NE.ll  f.fl  TO  9? 

WRITH3,93) 

93  FORMAT (  ,17X,7(  •  1 ' ,RX) I 

VRtTM3,97) 

97  FORMAT ! ,RX,73( 1 ) 

9?  IF  (  K3I  K5  l.L  T .  J1  )  GO  10  H7 
99  K*  =  K  1 ( K 9 1 

WRITF!  3,9ft) (AA( I  1, 1  =  1, KR1 ,AA(  109) 

K9=K9+l 

!F  (  K3I  K5  l.GE.  J 1  )  GO  TO  99 
«7  CONTINUE 
WRITE!  3,43) 

WRITE! 3  ,97) 

WHITH3.9M)  (  INTI  J).J  =  1.5) 

9«  FORMAT! 7X , 3 ( 1 3 . 1 5X 1 , 1 4 , 1 5X , I  3/ , IX ) 

WRITF(3,177) 

177  FOKMATI  34X , 'PLOT  OF  GO  VERM'S  THFTA'l 
IFIM.FG.l)  WR!TF|3,17A) 

17m  FORMAT!  '*•  > 

IMm.FO.?)  wilt  TF(  3,179) 

179  FORMAT  I  ' ♦ • ,4?X , > 

171  Cool  t imi if 
170  WRITE (3,104) 
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1.191 

-128.5 

1.219 

1.104 

27.4 

95.0 

1.217 

1.103 

-159.0 

1.171 

1.082 

-9.5 

100.0 

0.850 

0.922 

167.7 

1.122 

1.059 

-46.5 

106.0 

0.490 

0.700 

126.9 

1.072 

1.035 

-83.4 

110.0 

0.312 

0,559 

71.1 

1.024 

1.012 

-119.8 

115.0 

0.417 

0.646 

11.4 

0.977 

0.988 

-155.5 

120.0 

0.  791 

0.889 

-31.6 

0.930 

0.965 

169.9 

125.0 

1.220 

1.149 

-62.6 

0.883 

0.940 

136.7 

130.0 

1.848 

1.360 

-87.7 

0.831 

0.912 

105.3 

135.0 

2.232 

1.494 

-109.1 

0.770 

0.878 

76.0 

140.0 

2.387 

1.545 

-127.9 

0.697 

0.835 

49.3 

145.0 

2.296 

1.515 

-144.4 

0.608 

9.780 

25.2 

150.0 

2.  COO 

1.414 

-158.8 

0.505 

0.71 1 

4.0 

155.0 

1.575 

1.255 

-171.1 

0.391 

0.626 

“14.1 

160.0 

1.104 

1.050 

171.7 

0.275 

0.524 

-29.1 

165.0 

0.661 

0.813 

170.  7 

0.167 

0.408 

-40.6 

170.0 

0.306 

0.553 

165.0 

0.078 

0.280 

-49.3 

175.0 

0.078 

0.280 

161.5 

0.020 

0.142 

-54.3 

180.0 

0.000 

0.000 

0.0 

0.000 

0.000 

0.0 

PLOT  Of-  (RFAL  J2)/MAX>REAL  J?>  VERSUS  LENGTH 
PLOT  OF  MHAG  J2  J  /  HA  X  ’  I  M  AG  J«*  VERSUS  LENGTH 
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100.000  4- 


10.000  -I- 


APPENDIX  D.  PROGRAM  AND  SAMPLE  INPUT-OUTPUT  DATA  FOR  SECTION  VII. 

//  (0034.EE  ♦  8,2  >,  'MAUTZ,JOE*  ,MSGLEVEL=1 

//  EXEC  FORTGCLG,PARM.FnRT='MAP» 

//FORT .SYSIN  00  4 

SUBROUTINE  PLANE) VVR.THR, NT) 

COMPLEX  VVRU)  ,A5,A6,U 

COMMON  U,R(42),2S(42),SV(42),CV(42),BK,NP,NN,T(80),TR(fl0) 
DIMENSION  BJ(126),THR(1),FK(20) 

KG=NP-1 
NM=KG/2— 1 
.  M2=NN+? 

A5=2. #3*14159 3 MU**(NN+1) 

NV=NM*4 
FK ( 1 ) =1 . 

DO  153  J=1,M2 
J1=J+1 

FK( J1)=FK< J)*J 
153  CONTINUE 

00  156  L«1,NT 
L1=(L-1)*NV 
CS  =  COS ( THR  (  L ) ) 

SN=S I N ( THR ( L ) ) 

BCS=BK*CS 
DO  302  J=1 ,KG 
X=R( J)#BK*SN 
J1=J 
1 1  =  NN 

I F ( 1 1 )  303,304,303 

304  11=11+1 
J1=J1+KG 

303  00  305  JJ=I1,M2 
I F ( X-l *E“5 )  1,1,2 

1  IF(JJ-l)  3,3,4 

3  B  J  (  J  1 )  =  1 . 

GO  TO  306 

4  BJ ( J 1 ) =0 . 

GO  TO  306 

2  RH=X / 2 • 

RH2=RH*RH 

RM3  =  RH#*( JJ-1  ) 

BJ(J1)=RH3/FK(JJ) 

SS  =  B J ( J 1 ) 

8  SST=SS*l.+-7 
00  155  K= 1 ,20 
SS=-SS*RH2/K/(K+JJ-1 ) 

BJ( J1)=BJ( Jll+SS 
I F ! ABS ( S S )  —  S ST  )  306,306,155 
155  CONTINUE 
STOP  155 

306  J1=J1+KG 

305  CONTINUE 
302  CONTINUE 

IF(NN)  307,308,307 
30ft  DO  309  J= 1 ,KG 
J1=J+2*KG 
B J ( J ) =-B J ( J 1 ) 

309  CONTINUE 

307  00  300  J=1,NM 
J1  =  J  +  L  1 

J2= J 1+NM 
J3= J2+NM 


J4= J3+NM 
VVR ( J1 ) =0. 

VVR(J2)=0. 

VVR ( J3 ) =0. 

VVR< J4)=0. 

00  301  1  =  1,4 
1 1  =  2# l J-l )  + 1 
I4=4*( J-l )+I 
12=1 1+KG 
1 3= ! 2+K6 

A4=(C0S!ZS( I1)»BCS)+U*SIN(ZS( 1 1 ) *BCS ) ) *A5 
BJ1=(BJ< 1 3 ) +BJ  < 1 1 ))#.5 
BJ2=(BJ(I3)-BJ(I1))*.5 

VVR(  J1)=VVR(J1)+A4*!CS*SV(  II  )*BJ2+SN*CV(1 1)*BJ(  1 2  )*U )  =<«T  {  14) 
VVR!  J2)=VVR(  J2)+A4#CS*BJ1*U*TR( 14) 

VVR! J3)=VVR< J3)-A4*SV<I1)*BJ1*U*T! 14) 

VVR( J4)=VVR! J4)+A4*BJ2*TR( 14) 

301  CONTINUE 
300  CONTINUE 
154  CONTINUE 
RETURN 
END 

SUBROUTINE  REORO(K1,K3,L  ) 

DIMENSION  K1(1),K3<1) 

DO  81  J  =  1  ,L 
K8=K3( J  ) 

K4= J 

DO  82  l  =J»  L 
I F ( K3 ( I ) -K8  )  82,82,84 
84  K8=K3( I  ) 

K4=  I 

82  CONTINUE 

K3(K4)*K3< J) 

K3 ( J ) =K8 
K8=K 1 ( K4 ) 

K1!K4)=K1< J) 

Kl ( J ) =K8 
81  CONTINUE 
K3 ( L+l ) =- 1 
RETURN 
END 

COMPLEX  A3, Y ( 1 400 ) , VVR < 5840 ) ,63( 40 ) ,64 ( 40 ) , E l ( 73 ) , E? ( 73 ) »U 
COMMON  U,R(42) ,ZS(42),SV(42),CV(42),BK,NP,NN,T(80),TR(R0) 
DIMENSION  RH(43),ZH(43),DH(42) , T J ( 20  ) , I  NT ( 5 ) , THR ( 73 ) 
DIMENSION  AA< 110) ,K1(73> ,K2!73) ,K3<73) ,K4(73) 

DATA  AA(1),AA(1Q9),AA(H0)/'  ',»X',*0'/ 

DO  107  1=1,107 
107  AA( I+l )=AA! I ) 

U=(0.,1.) 

ET  A=37  4,707 
PI=3. 141593 
PR=180»/PI 
REWIND  4 

50  READ! 1,51, END=52)  KK,NP,NT,BK 

51  FORMAT(3I3,E14.7) 

READ! 1,53) (RH( I ) ,1=1, NP) 

READ! 1,53) (ZH( I ) ,I=1,NP) 

53  FORMAT! 10F8, 4) 

74  WR I TE ( 3 ,54 )  KK ,NP »NT ,BK 
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54  FORMAT  (  IX//  •  KK=',I3,'  NP=',I3,'  NT=»,!3,»  ftK=',E14.7) 

WR I TE ( 3  *55 ) 

55  FORMAT (IX/*  RH* ) 

WRITE(3,46) (RH(I) ,I=1,NP) 

46  FORMAT ( IX, 10F8.4) 

WR I TE ( 3  *56 ) 

56  FORMAT( IX/'  ZH» > 

WR I TE ( 3,46  )  (  ZH  (  I ) ,!  =  1,NP) 

KL=  1 

126  1F( (RH(l)-RH(NP) ).NE.O..OR.(ZH(l)-ZH(NP)).NE.O.)  00  TO  5« 
KL  =  0 

RH ( NP  +  1 ) =RH { 2 ) 

ZH( NP+ 1 )  =ZH( 2 ) 

RH ( NP+? ) =RH ( 3 ) 

ZH ( NP+2 ) =ZH( 3) 

NP=NP+2 

58  00  57  1=2, NP 
12=1-1 

RR  1=RH(  I  )-RH< 12) 

RR2=ZH( 1 )  -ZH  { 12) 

DH( 12)=S0RT(RR1»RR1+RR2»RR2) 

ZSU2)  =  .5*(ZH(1)+ZH(  12)1 
R  (  I2)=.5*(RH( I ) +RH I  12)) 

SV( I2)=RR1/DH( 12) 

CV(  I2)=RR2/l)H(  12) 

57  CONTINUE 

l)T  =  Pl  /  (NT-1  ) 

00  1  J=1 ,NT 
THR(J)=DT*( J-l) 

1  CONTINUE 
BK2=BK*BK 

Pl=S0RT(BK248K2*ETA*ETA/4./Pl»*3) 

NM=(NP-3>/2 
NM4=NM*4 
NM2=NM*2 
NZ=NN2*NM2 
DO  74  1 , NM 

J2=2*  (  J-l )  ■►  1 
J3  =  J2+  1 
J4=J3+1 
J5= J4+ 1 
J6=4*(  J-D  +  l 
J7= J6+1 
J8=J7+1 
J9  =  J8+1 

DELI =DH ( J2) +0H ( J3 ) 

OEL2=OH( J4)+0H{ J5) 

TC J6)=0H( J2)*0H( J2)/2./0Ell 
T  (  J7 ) =0H ( J3 ) *( DH( J2 ) *DH( J3 ) /2» ) /DEL  1 
T( J8)=0H( J4)*(DH( J5)+DH( J4)/2.)/0EL2 
T(  J9)=0H( J5)*OH(J5)/2./OEL2 

74  CONTINUE 

00  75  J= 1 , NM4 
TR ( J ) =T ( J ) 

75  CONTINUE 

115  IF(Kl.EO.O)  GO  TO  78 
I F ( RH ( 1 ) )  77,23,77 
77  DEll=OH( 1 ) +DH ( 2 ) 

TRC1)=0H(1 )*<1.+(OH(2)+OH(1 )/2. )/DELl ) 
TR(2)=OH(2)*(l.+OH(2)/2./DELl) 
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23  !  F ( RH( NP ) )  79,78,79 
79  J 1= ( NM-1 ) *4+3 
J2=J1+1 

DEL2=DH(NP-2)+DH'NP-l) , 

TR1 Jl)=0H(NP-2)*( l«+OHlNP-2)/2*/OEL2) 

TR(  J2)=0H<NP-1)*(  l.+(DH(NP-2)+DH(NP-l)/2.)/0EL2) 

78  SS=0. 

00  7  1=1, NM 
11=2*1 I-1J+1 
12=11+1 

SS=SS+0H{ 1 1)+0H( 12) 

TJl I )=SS 

7  CONTINUE 
DEL=TJ(NM) 

IF1KL.NE.0)  DEL=DEL+DH<NP-2)+DH<NP-l) 

OEL=OEL / 10. 

OH  8  J= 1 *NM 
TJl J)=TJ( Jl/OEL 

8  CONTINUE 

DO  42  J=1,NT 
Ell J)=0. 

E2 1 J ) =0 . 

42  CONTINUE 

00  40  M=1,KK 
NN=M-1 

CALL  PLANE! VVR,THR,NT) 

127  READ! 6 ) 1 Y< I ) , 1=1, NZ ) 

U2  F0RMAT?/3x!»0'  ,4X,  •  SIG  00-,2X,-MAG  SOO*  ,1X,  -SIG  00-,?X,-MAG  800- ) 

113  FORMAT ( ' +  -• *8X, * — *,7X,» — 5X 7X, •//• ) 

1)0  41  L  =  1  ,NT 
L1=(L-1 )*NM4 
Of)  2  J  =  1,NM2 
E3 1 J ) =0* 

E4 1 J ) =0 . 
on  3  1=1, NM 
Jl=J+< 1-1 }*NM2 
J2=Jl+NM*NM2 
1 1  =  1 +L  1 
12= I l+NM 
13=11 +NM2 
1 4= l 3+NM 

631  J)=E3(  J)+Y(  J1)*VVR! Il)-Y(  J2)*WRU2) 

E4( J)=E4< J)-Y( J1)*VVR1 I31+Y1 J2>*VVR( 14) 

3  CONTINUE 
2  CONTINUE 

on  43  J=l»NM2 
J 1 = J  +  L  1 
J2= J 1 +NM2 

E1!L)=E1(L)+E3( J)OVVR(Jl) 

E2(L)=E2(L)+64< J)*VVR(J2) 

93  CONTINUE 

JFtNN)  44,45,44 
45  El  1  L  )  =  ,5*t !  <U 
E2(L)=.5*E21L> 

44  X1  =  THR(L  )*P'< 

X3=Pl*CAKSiei(D  ) 

X2=X3*X3 


X5=Pl*CAHS(h2(L)  ) 

X4=Xb*X5 

WKIThr3.ll])  XI  ,X2,X3,X4,Xb 
111  H1*«AT(1X,h5.1.4HH.3) 

41  Cl'NT  I  t'lllh 
40  CONI  1  rytM- 

W«  1 T  b  C  i,t<)6) 

w«  i  n.  { 3, 1 1?) 

WK!  II- (3, 113) 
on  an  j=i,nt 

K1 (J)=TMK ( J)#7?./PT*8.5 
K2( J)=K1 (J) 

K3t  J)=20.*AL0G10IPl*CAH<i(Kl<  J)  ) )+?0.5 
K4 ( J ) =20.4ALUG10( PI *CABS ( 6? ( J ) ) ) +20. 5 
X1=THKU)*PR 
X3  =  P1*CAHS(E1( J  )  ) 

X2  =  X3<=X3 

Xb  =  Pl*CAnSO:2<  J  )  ) 

X4=X5*Xb 

WPITK3.111)  XI  ,X2,X3.X4,Xb 
flu  COM  iNUh 

14  CALL  KhOKlMKl.K3.NT  ) 
i*>  CALL  KHIKIMK2.K4.NT  ) 

DO  10b  J  =  1 ,  b 
INI ( J ) = ( J-l ) *45 
10b  CUNT  I Ntlt 
XI  =  moo. 

KS=1 
K6=  1 

WK  I  TH  3,  106) 

106  MIKMAT  (  •  1  •  ) 

01)  «7  J  =  l  ♦  *5 1 

Jl=Sl-J 

WK1TK3.HH) 

AA  FMRi»AT(9X.*  I  •  ,71  X,  •  I  •  ) 

!H  t J-l )/ 10*10- (J-l ) )9?,40,«? 

40  WK1TH3.9})  XI 

41  hOHMAU  ,h7.2.  *  — 1 , 64X  ,  1  —  •) 
xi=xi/ln. 

1H  J.Nh.l  )  fill  TO  42 
WK  I  T  M  3,93 ) 

93  H1KMAT ( , 17X.7I • l * ,KX) ) 

W«  I  rt ( 3,4  7) 

97  H1KMA1  ('■*•'  ,  H  X  ,  7  3  (  J  > 

42  IHK3IKS)  .LT.J1  )  GO  TO  94 
4 5  KH  =  K  1  ( Kb  ) 

WK  ITfc  f  3.4K 1  ( A A l  '  )  ,  I  si  .KH  )  ,AA( 109) 

4 A  M1MMAT  (  '  ♦»  ,  1  1 0 A  l  1 
K  b  =  K  b ♦  1 

IP  (K3f  Kb  )  .Gb.Jl  )  0(1  TO  95 
44  IMK4IK6)  .LT.J1  )  0(1  TO  H7 
96  KR-K2(K6) 

WKITfc 13,4 A | (AA(l) ,1  =  1 ,KH) ,AA(1  10) 
K6=K6+ l 

IF (K4{K6)  .Ot  .J1  )  OH  TO  96 
47  CONIINUh 
WKITH  3,93) 

WK I Th ( 3,97  1 

WKI  IH  3,9H  )  (  1NT  I  J)  ,  J=1  ,b) 

9 A  HIKMAT(7X,3(I3,lbX),I4,15X,!3/.lX) 
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WRITE (3.99 ) 

99  FORMAT ( 16X, • X  X  X  PLOT  OF  SIGMA  00  OVER  LAMBDA  SQUARED  VERSUS  THE 
IT  A  •  ) 

WR I Te  <  3, 101 ) 

101  FORMAT!  •♦S36Xt»  — •  ) 

WR I TE ( 3» 100 ) 

100  FORMAT ( 16X, *0  0  0  PLOT  OF  SIGMA  00  OVER  LAMBDA  SQUARED  VERSUS  THE 
IT  A  *  ) 

WR 1 TE ( 3* 102  I 

102  FORMAT ( '+• ,36X, •//•  ) 

WR l TE ( 3 , 106 ) 

GU  TO  50 
52  STOP 
END 

/* 

//G0.FT06P001  OD  USNAM6=EE0034. REV1 , 01 SP  =  0LD,UNITn;>314,  X 

//  VOLUME=SER=SUOOOA, DCHs (RECPM=V*BLKSUE=1R00,LRECL= 1796) 

//GO. SYS  IN  DD  * 

00509 1037  0.4659995E+00 


0.0 

0.0868 

0.1736 

0.2605 

0.3473 

0,4341 

0.5209 

0.6078 

0.6946 

0.7814 

0.8682 

0.9551 

1.0419 

1.1287 

1.2155 

1.3024 

1.3892 

1.4760 

1.5628 

1.6497 

1.7365 

1.8233 

1.9101 

1.9970 

2.0838 

2.1706 

2.2574 

2.3442 

2.4311 

2.5179 

2.6047 

-0.0000 

2.6837 

2.6863 

2.5969 

2.4184 

2.1570 

1.8216 

1.4238 

0.9772 

0.4971 

0.0 

0.4924 

0.9848 

1.4772 

1 .9696 

2.4620 

2.9544 

3.4468 

3.9392 

4.4316 

4,9240 

6.4164 

5.9088 

6.4013 

6.8937 

7.3861 

7.8785 

8.3709 

8,8633 

9.3557 

4.8481 

10.3405 

10.8329 

11.3253 

11.8177 

12.3101 

12.8025 

13.2949 

13.7873 

14.2797 

14.7721 

18.2260 

15.2657 

15.7650 

16.2562 

16.7225 

17.1478 

17.5177 

17.8195 

18.0427 

18.1790 
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KK=  5  NP 

=  41  IS  T  = 

37  BK  = 

9. 465944%  jO 

8H 

J  •  V 

C.7136d 

C. 1736 

C.  2  605 

C„  >4 73 

L  .  4  34  1 

.52  9 

'./.r  ?n 

**  ,694/. 

’  .7814 

C.8652 

C.9551 

1.3414 

t. 1287 

i.2-55 

1 . 3024 

1.3892 

1  .^76 4 

1 .46,. '9 

!  .  /  44  7 

1.  73t’6 

1.8233 

1.9111 

1.9973 

2*Cf38 

?.  17  6 

?.?574 

2.  14<.  ? 

2.4  31  1 

,'.4  179 

2.604? 

-0.0 

2.6837 

2.4863 

2.5964 

2*^i 84 

2. 157C 

1.8216 

1.4238 

1  .4  7  72 

• .49/1 

ZH 

o.c 

C .4024 

7 .0848 

1.4  7  72 

1.9646 

2.462C 

2.9544 

3.4468 

3.9342 

4,4316 

4.0240 

5.4164 

6. 4088 

6.4  31  J 

6.  8437 

7. 3861 

7.8  785 

0  .  3  7*'  9 

i.'ltill 

1.  144  7 

9. 84u 1 

1C.  14(5 

1C .8329 

11. 3253 

11.8177 

1?.  ii  .1 

12.8025 

13.2949 

13.7873 

14.2747 

14.7  7?1 

18.226C 

1 5. 2667 

15. 7640 

16.266? 

16. 72?5 

17. 1478 

1  7.41  77 

17.K195 

18.  4>7 

ll>.  !  748 

8 

SIG  68 

MAG  S8P 

S1G  1? 

MAG  S20 

u.:' 

O.C 

■ 

0.0 

i.O 

5.r* 

o.rc  1 

).i  25 

C  .07>t 

G.0C1 

10. r 

o.'h: 

J.K2 

2.000 

6. 004 

15.. 

0.151 

3.225 

C.0C1 

C.O  1> 

2C  . . 

u.  152 

i,  .  389 

C.CoO 

C.017 

25. C 

0.337 

0.681 

0.001 

C.027 

3t.'. 

6 . 6C  3 

0.776 

C.0C1 

0.038 

35.  v 

u  .  84  ■ 

c .  00  3 

52 

4<  .' 

l.C4<, 

1.7  44 

:.C05 

0.068 

4  5.'. 

1.  107 

1 5? 

C  .007 

0.087, 

61  .  , 

0.9U 

c .  0 1 1 

9.105 

55. ; 

U.  6 1C. 

).  781 

0.015 

0.124 

69. t 

0.352 

'*.64  J 

0?1 

144 

65. C 

0.  182 

..427 

0.7130 

174 

7  O.C 

0  .f  53 

' .  ?  3v. 

-  .052 

0.228 

75. 0 

P.C06 

G  .7  80 

0.095 

0.3C8 

ao.c 

0.253 

0.5C  3 

0,158 

0.397 

h5. r 

0.891 

•3.944 

7!.  221 

0.47'' 

9G.C 

1,605 

1.26  7 

0.254 

G.5C4 

95.6 

1,871 

1.368 

0.238 

C.4HH 

1  oc . » 

1.499 

i,?24 

C  .  1  8 1 

G.  426 

105,0 

O.PU 

0.40? 

C.  112 

0.335 

1  10. c 

r .  266 

0.51/. 

0.056 

0.237 

1  1 5 .  C 

C.C33 

0.  181 

O.C23 

C.  153 

12C  .0 

C.CK 

: .  i  c  1 

0.  009 

C.094 

125.C 

C.C43 

0.208 

0.004 

0.062 

1  3C.Q 

C.062 

0.250 

0.002 

0,047 

1  35.0 

C.C65 

0.254 

0.031 

0.037 

140. C 

0.C59 

C  .  ?4  3 

0.001 

0.028 

145.0 

0.C44 

0.22? 

G  .000 

0.0271 

1  50. C 

o.C  36 

o.iaa 

0.000 

0,014 

15  5.4 

0.C21 

0.146 

0.000 

C.010 

160.0 

0.010 

0.101 

o.oco 

0.006 

164.0 

0.C04 

rt.C  60 

C.000 

0.003 

1  70.0 

o.cot 

0.C27 

0.000 

0.001 

175.0 

0.000 

0.CC7 

0.000 

0.000 

180.0 

o.cco 

0  .  c  00 

c.o 

0.0 

8 

sic  ee 
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